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SOME MINERALOGICAL EQUILIBRIA 
IN THE SYSTEM K,0—AI,0;—Si0,—H,0 
J. JULIAN HEMLEY* 


Department of Geology, University of California, Berkeley, California 
ABSTRACT. An experimental investigation was made at elevated temperatures and 
pressures of the hydrolysis of K-feldspar to mica + silica and of mica to kaolinite, Values 
of the equilibrium quotient, mxei/mne:, were determined at temperatures ranging from 
200° to 550°C. The solution pressure for most of the runs was 15000 psi although deter- 
minations were made at both higher and lower pressures, The reactions are exothermic 
and the equilibria show a marked shift to higher acidities with increasing temperature. 
The equilibrium curves for the two reactions are essentially parallel, Above about 350°C 
the assemblage mica-kaolinite changes to mica-pyrophyllite-boehmite, and at still higher 
temperatures mica-pyrophyllite-andalusite is formed. At 400°C and 15000 psi the equilib- 
rium quotient is 10°" for the assemblage K-feldspar-mica-silica, and 10°* for the assem- 
blage mica-pyrophyllite-boehmite. The experimental equilibrium quotients probably do not 
differ from the thermodynamic equilibrium constants of the reactions by more than a 
factor of 5 up to temperatures of 350°C or higher. 

The experimentation shows that the most important controls on the fields of stability 
of these several minerals are the K*/H* activity ratio and the temperature. At a given 
temperature and with increasing K*/H* ratio the fields of kaolinite, mica and K-feldspar 
are successively traversed, Similarly, at constant K*/H* ratio with increasing temperature 
this same sequence is observed. The free energy drive of hydrolysis decreases with in- 
creasing temperature. Thus the alteration potential of hydrothermal solution would increase 
with decreasing temperature provided the rate of migration were sufhiciently rapid in 
comparison with the rate of reaction with the wall rock. Experiments at solution pressures 
of 5000 and 35000 psi show that the effect of solution pressure on these equilibria is rela- 
tively small. An increase in pressure increases the extent of hydrolysis, whereas a decrease 
in pressure shifts the equilibria in the opposite direction. 

Phe experimental results have particular application to hydrothermal alteration and 
the genetic interpretation of various alteration mineral patterns. 


INTRODUCTION 
Detailed examination of many occurrences of rock alteration in ore de- 


posits over the last 20 years has emphasized the need for a greater understand- 
ing of the chemistry of rock alteration processes. Advances in the study of 


sheet-structure silicates during this period have made possible the identification 


and structural characterization of minerals composing rock alteration suites, A 
number of geological investigations of rock alteration have utilized and con- 
tributed to this fund of information. These studies have revealed systematic 
features, as well as complexity, in the mineralogy and field patterns of mineral 
distribution. They have thus brought into focus a great lack of fundamental 
chemical information on alteration processes and the genetic significance of 
various alteration patterns. 

The term rock alteration as used here has the general meaning associated 
with it in the geological literature regarding ore deposits, It is characterized 
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chiefly by the development of primarily layer-structure, hydrated minerals 
from pre-existing rock forming minerals, notably feldspars. 
The present study considers equilibria in the system K.O—AI,0,—SiO 
H.O. Consideration of these components involves, among other minerals, the 
stability relations of sericite, the fine-grained potassium-aluminum mica. 
Sericite is one of the most characteristic alteration minerals associated with 
hydrothermal sulfide ore deposits. A better understanding of the chemistry of 
sericite will have broad application to alteration problems in general. More- 
over, sericite, as: compared to other typical alteration minerals, is most often 
closely associated spatially and temporally with the ore itself. Information on 
the chemistry of sericite will therefore relate closely to the problem of the 
nature of the processes by which ore minerals are deposited. 


SOME FIELD RELATIONS 

This investigation was the outgrowth of studies of rock alteration in hy- 
drothermal vein deposits in granodiorite-type rocks, and alteration in the 
porphyry coppers. The geological aspects of the study have been summarized 
in another paper (Meyer and Hemley, 1957) and are only briefly noted here. 

Perhaps the most unambiguous geometrical and temporal relationships 
between alteration mineral assemblages are observed in certain vein deposits 
of the type mentioned above, notably those of Butte, Montana, Enveloping the 
copper, zinc and manganese veins of Butte an alteration envelope made up of 
two major units is observed. An innermost sericitic envelope is bordered sharp- 
ly by an outer argillic envelope which in turn grades into fresh granite, The 


clay mineral content of the argillic envelope consists of a high percentage of 


kaolin nearer the vein with montmorillonite dominant farther out, The de- 
velopment of the clay zone is the result of the preferential decomposition of 
plagioclase. K-feldspar and quartz are only slightly affected. In the sericite 
zone quartz and white mica are the dominant mineral constituents, Clays are 
absent. K-feldspar in this zone in some environments, such as in the upper 
levels of a large portion of the Butte district, is completely sericitized, In othet 
instances, notably in some alteration types in porphyry copper deposits, the 
association of K-feldspar with sericite, as well as with clays, is commonly 
obse rved, 

These brief summary statements should not be considered an adequate 
review of all aspects of the rather variable alteration picture under considera- 
tion. They do, however, serve to delineate a number of important chemical 
problems. May the assemblages mica-K-feldspar, mica-kaolin, and K-feldspar- 
kaolin, each in association with quartz. represent stable associations? If so, 
what are the parameters which define their stability? If, on the other hand. 
any of these assemblages are metastable, what kinetic factors are involved in 


their formation? Investigations of the reactions which give rise to these 
various mineral associations are requisite to an understanding of the conditions 
of their formation and to the definition of their respective stability fields, Of 


immediate interest are the hydrolysis equilibria of K-feldspar and K-mica, 


PREVIOUS WORK 


Exploratory experiments on the decomposition by feldspar in aqueous 


in the System K,O—AL,O,—SiO,—H,0 243 


solutions were carried out by Gruner (1944). The studies involved the treat- 
ment of feldspar in a large excess of 0.1N HCl at elevated temperatures and 
pressures, Some experiments were also conducted using chemical mixtures of 


equivalent feldspar or mica composition, The temperatures of experimentation 
were 300, 340, 360, and 400°C, Pressures in the bombs were determined by 
degree of filling. 

In using a large excess of acid, the experiments were designed so that the 
pH at the end of the run was approximately the same as at the beginning. In 
all instances in which K or Na feldspar alone were treated in the acid solutions 
pyrophyllite was produced. If aluminum hydroxide and KCl were added to 
the system, mica was formed at higher temperatures and kaolinite or boehmite 
at lower temperatures. The transition occurred between 340° and 360°C, If 
KCI alone were added in the treatment of feldspar no alteration was observed. 
It was concluded that under these conditions K-feldspar persists metastably. 
The results of Gruner’s work indicate a course of alteration under conditions 
of constant high acidity at certain conditions of temperature and pressure, but 
they have little import from the standpoint of alteration under equilibrium or 
near equilibrium conditions, They therefore have limited application to the 
definition of K-feldspar or mica stability fields. 

Studies of mineral alteration in a leaching environment under conditions 
of moderate temperature, pressure and acidity were conducted by Norton 
(1939). These experiments, like those just described, were designed as ex- 
ploratory tests on mineral alteration rather than as chemical equilibrium 
studies. The experimental technique involved continuous leaching of finely 
ground feldspar or other silicate mineral by CO,-bearing aqueous solutions. 
Sericite and quartz were reported to have formed from K-feldspar. Mineral 
identifications as they were reported, however, were not completely satisfactory. 

Morey and Chen (1955) conducted similar studies using pure water, The 
relative rates of solution of mineral components were determined and the resi- 
dual mineral substances identified. 

There has been considerable experimental work on synthesis of clay min- 
erals and mica under hydrothermal conditions. Much of this work, such as that 
of Noll (1936), has been adequately summarized in a comprehensive review 
by Morey and Ingerson (1937). An important investigation of mineral stability 
relations in the system K,O—AI,0,—SiO.—H.O was that of Yoder and 
Eugster (1955) on synthetic and natural muscovite. The studies of Yoder and 
Kugster, however, were concerned primarily with the thermal decomposition 
reaction of muscovite. Equilibria such as those involved in the present study 
were not considered. 

On the basis of the work of Gruner and others some investigators have 
attempted to deduce fields of formation for various alteration minerals under 
hydrothermal conditions (Folk, 1947; Stringham, 1952). Results of the pre- 
sent work indicate these proposed relations to be substantially incorrect, 


EXPERIMENTATION 
The foregoing discussions have indicated the need for investigation of the 
hydrolysis reactions of K-feldspar and K-mica. In view of the lack of appropri- 
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ate experimental data three fundamental points must be considered: (1) What 
are the correct net reactions involved? (2) What are the equilibrium constants 
of the reactions? (3) What is the significance of kinetic factors, reaction rates 
and mechanisms. in the occurrence of these reactions? The first two of these 
questions are logically of more immediate or fundamental importance and it is 
with these problems that the present study is largely concerned. 


Procedure 
Initial « xperimentation established the occurrence and reversibility of the 


reactions, 


K-feldspar mica 


>KAISi,O. + H + 3 SiO, + K 


mica 


KAI. AISi,O,,(OH). + H* = -yALSi,0,(OH), + K 


These experiments were conducted using dilute HCI-KCI solutions, For ex- 
ample, the spontaneous ilteration of K-feldspar to mica plus cristobalite and 
of mica to kaolinite can be accomplished in solutions of 0.15m HCI-0.5m KCL." 
These reactions are suggested by geological and chemical considerations, but 
the decomposition process for K-feldspar previously had not been experimental- 
ly realized. It is also conceivable that the decomposition of K-feldspar to mica 
would not represent an equilibrium process under hydrothermal conditions ap- 
propriate to rock alteration. The necessity for demonstrating a given net reac- 
tion can hardly be overemphasized. 

Following such studies of the nature of the reactions, equilibrium points 
were fixed more precisely by successive equilibration runs. Such experiments 
could conceivably be conducted in a number of ways. but with the equipment 
available to this project the most appropriate method appeared to be the 
equilibration ol an excess ol the mineral assemblage under consideration with 
an aqueous phase, followed by a measurement of the potassium and acid con- 
centration in that phase. Equilibrium was approached from hoth directions. 
that is, by leaching and by fixation of potassium, For example. in K-fixation 
experiments involving the K-feldspar-mica-silica assemblage. K-feldspar was 
produced from mica plus silica in the presence of aqueous potassium at an 
acid concentration initially less than the equilibrium acid concentration, K- 
leaching experiments corre spond d to the reverse process. 

The experiments were conducted using sealed platinum tubes, One hun- 
dred to two hundred milligrams of the appropriate mineral assemblage (< 200 
mesh) were used together with about a milliliter of solution of appropriate 
normality in HCl and KCI. The tube was placed in a cold-seal rod bomb under 
external water pressure and the bomb heated to temperature, The use of such 
vessels in hydrothermal experimentation has been described by Tuttle (1949) 

The alteration of K-feldspar to mica silica (cristobalite) in dilute KCI solutions in- 
dicates Gruner’s conclusions that feldspar remains unaltered under such conditions to be 
mistaken. On the basis of some preliminary results of the present study it is suggested that 


there was alteration of the finer grained material in Gruner’s charges, but the product 
was amorphous material not evidenced in his x-ray diffraction patterns. 
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and others.* The total confining pressure for most of the runs was 15000 psi. 
Temperature was controlled within +5°C. 

At the end of the run, the bomb was quenched rapidly at essentially 
constant pressure. The quench consisted of immersion in a mercury bath, con- 
tained in a large reservoir of ice-water mush. The bomb reached a low internal 
temperature, probably less than 50°C in 1-144 minutes, and then about 2 
minutes was required for extraction of the sample and measurement of its pH. 
Since the solutions contained KCl in large excess, the concentration of this 
material was known. The concentration could of course be calculated in those 
instances where there was significant change in KCl and HCl concentration in 
the experiment, The true HCl concentration of the experimental samples was 
readily obtained by calibration of the response of the pH meter with solutions 
of known HCI-KCI normality. A Beckman one-drop glass electrode was used. 
The reproducibility of the readings was +0.03 pH units. 

Both synthetic and natural minerals were used in the experimentation. 
The same equilibrium results were obtained with each, equilibrium being ap- 
proached from either direction, Synthetic materials yielded slightly faster re- 
action rates and were used in a majority of the runs, The synthetic minerals 
were made hydrothermally from gels of appropriate composition, Reaction 
time for the runs was from one to three weeks. At lower temperatures experi- 
ments of several weeks duration were the rule. Approximately one hundred 
runs were involved in the experimentation. 


EXPERIMENTAL EQUILIBRIUM CURVES 


Mica — Pyrophyllite — Andalusit 
K-feldspar— Mica —Silica 
Mica — Pyrophyllite —Boehmite 


Mico — Kaolinite 


S 
K-feldspor 
Mica 
Py. - \ 
Boeh ~ 


\ \ 
WA 
@ fixation of K 


Oo leaching of K 


15000 ps: total pressure 


10" io? 
MHCI 


Fig. 1. The experimental equilibrium curves. K-fixation designates the approach to 
equilibrium from the right on the diagram, K-leaching designates the reverse process. 


* The writer is grateful for two bombs loaned by George Kennedy as models for the ap- 
paratus, 
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Experimental Results 


Stability Fields.—Results of the experimental work are summarized in the 
equilibrium curves of figure 1. The pressure of the solutions was 15000 psi. 
The curves delineate the fields of stability of K-feldspar, mica and kaolinite. 
They are expressions of the critical KCl/HCI ratio for the transition from one 
field of mineral stability to another at a given temperature, The KCI/HCI 
ratio is an approximate expression of the equilibrium constant of the reactions, 
and shows a shift of the equilibrium to more acid conditions at higher tem- 
peratures. The KCl/HCI ratio, the temperature and the pressure are therefore 
the parameters which define stability relations among these mineral assem- 
blages. 

The variation of the equilibrium quotient with temperature for the two 
mineral assemblages, K-feldspar-mica-silica and mica-kaolinite, are given in 
figure 2. The curves are essentially parallel and separated by an interval of 
acid concentration corresponding to about 1.5 pH units. The significance of 
the dashed curve is discussed on a later page. At a given KCI/HCI ratio, then, 
the stability fields of kaolinite, mica and K-feldspar are successively traversed 
with increasing temperature. Similarly, at a given temperature the three fields 
are successively traversed with increasing KC]/HCI ratio. 

At higher temperatures the curve corresponding to the mica-kaolinite 
equilibrium gives way to different relationships. At 350-355°C mica de- 


Fig. 2. Plot of the experimental equilibrium quotients versus temperature. The dashed 
curve is ussed on i iater 
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composes to pyrophyllite plus boehmite instead of kaolinite because of the in- 
stability of kaolinite above this temperature. 


mica  pyroph. ,  boeh. 


KAI. AISi,O,.(OH). + H+ = + + K+ 


kaolin pyroph. boeh. 


or Al,Si,0; (OH), ALSi,O, (OH). + + 


The point at 352°C on the 0.65m curve shows only minor decomposition of 
kaolinite whereas on the other curves kaolinite is replaced by pyrophyllite + 
boehmite. Diaspore rather than boehmite may be the stable alumina hydrate 
under these conditions (Roy and Osborn, 1954), but this would have very 
little effect on the energy relations. 

At 400°C the assemblage mica-pyrophyllite-boehmite was observed and its 
equilibrium quotient determined. At somewhat higher temperature the as- 
semblage mica-pyrophyllite-andalusite is developed. 


mica pyroph. andal. 


KAI. + + Al.SiO; + H,O + Kt 


The temperature of formation of this assemblage, as well as the dehydration 
temperature or upper stability limit of the alumina hydrate, are among con- 
siderations that could not be investigated fully in the present study because of 
limitations of time. The tendency of the monohydrate to persist above its tem- 
perature limit of stability has been noted by other investigators (Roy and 
Roy, 1955). The temperature of formation of the mica-pyrophyllite-andalusite 
assemblage is at least as low at 450°C. Evidence of the presence of andalusite 
in some experiments at 420°C suggests that andalusite will form at this tem- 
perature under the conditions of the present experimentation if sufficient time 
is allowed for crystallization, The development of andalusite agrees with the 
work of Clark, Robertson and Birch (1957) on the stability fields of Al,SiO; 
polymorphs. The x-ray powder diffraction pattern of the andalusite phase is 
nearly identical with that of andalusite from Fresno County, California. 

As discussed in an earlier paragraph the equilibrium points illustrated in 
the curves in figure 1 were determined by the degree of approach of potassium 
leaching and potassium fixation experiments. In K-leaching the approach to 
equilibrium was from the left on the diagram; for K-fixation it was from the 
right. For some of the equilibrium values a number of successive experiments 
were required before the two approaches closed on one another in the time 
allowed for reaction. In other cases, such results were obtained in one or two 
runs. The reason for a lag in approach to equilibrum in many instances is 
presumably grain-coating effects. When the position of the curve becomes ap- 
proximately known, experiments can of course be run with certain attainment 
of equilibrum. The values shown in figure 1 represent nearly all of the equilib- 
rium points. Experiments to determine the reactivity of dilute HCl solutions 
with platinum indicated the reaction to be negligible. An experimental un- 
certainty in this experimentation of 0.15-0.20 pH units is indicated by typical 
equilibrium values plotted on the curves. 
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The x-ray diffraction results for a number of the runs are shown in figure 
3. Powder diffraction patterns were of course the principal means of mineral 
identification because of the very fine grained character of the reaction prod- 
ucts. The muscovite polytype formed in the experimentation was generally 
1Md, the l-layer disordered structure. The 2M, polytype was indicated in a 
few runs. 

Slow crystallization of silica at low temperatures.—In the K-feldspar-mica- 
silica assemblage at temperatures below 400°C, considerable trouble was en- 
countered because of the tendency of silica to remain uncrystallized. The result 
of this was a higher activity of silica in the solutions than that present at 
equilibrium, with consequent stabilization of K-feldspar at acidities much 
higher than the equilibrium value. This supersaturation with respect to silica 
was markedly illustrated by experiments utilizing powdered hydrated silica 
(reagent silicic acid). At 400°C and above, this material readily crystallized 
to cristobalite or silica-K, but at 300°C and below, there was no crystallization 
of a silica mineral. Instead, the reacted solutions, after extraction from the 
tubes and standing for several days, congealed to a gel. 

The problem presented by the slow crystallization of silica required that 
the equilibrium point in these instances be established primarily by K-fixation 
using quartz, Even synthetic quartz, however, because of its extremely fine 
crystallinity and therefore greater solubility was possibly a source of slight 
error. The points on the K-feldspar-mica-silica curve at 200°C represent syn- 
thetic quartz, on the one hand, and natural quartz sized to >10 microns on 
the other. The sized material was slow in reacting and the result probably does 
not express complete approach to equilibrium. If the natural material had not 
been sized it would probably have given a result falling essentially on the 
curve as drawn, At higher temperatures where silica crystallized readily, the 
use of either cristobalite or quartz in the experiment yielded the same equilib- 
rium results within experimental error, This would be expected from the small 
difference in free energy between these two phases. 

Hydrolysis on Quenching.—Referring to the experimental curves it is 
seen that the equilibrium acid concentration in the several reactions is in- 
creased by higher temperature, as well as by greater salt concentration, These 
higher concentrations of salt and acid would be advantageous in minimizing 
errors that might result from back hydrolysis during quenching. Back hy- 
drolysis reactions would probably have to be largely of surface character in 
order to occur to a significant extent during the short time interval allowed 
in quenching. The amount of reaction that could occur would therefore be 
limited, and under more acid conditions the resulting error would be less 
significant. 


If the equilibrium acidity corresponded to a high pH value, say 4 or 5, 
only a negligible amount of reaction, essentially 0.0001 meq/ml, would be suf- 
ficient to change the pH by a whole unit. On the other hand, in order to 
change the pH by a whole unit in a solution with an acid concentration cor- 
responding to a pH of 1, the amount of reaction required would be about 0.1 
meq/ml. This is a very large amount of reaction, of the order of the total 
amount of mineral present in the charge. Also, it cannot be argued that higher 
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acidity would necessarily increase the rate of hydrolysis since increased 
acidity implies higher salt concentration which would suppress hydrolysis. 
These statements are born out by observation, In some of the initial experi- 
mentation. a slow drift in pH to higher values aftet quenching Was sometimes 
observed, but only in those solutions with pH values greater than 4 or 5. 

The rel tionship between free energy and rate of reaction is also of sig- 

nificance. From the slope of the curves it is apparent that the greatest free 

hydrolysis would exist at low temperatures, where rates of 

ould be minimum, At high temperature, where rates are high, the 

tendency hvdrolysis diminishes. It is therefore possible that hydrolysis 

during quen hin would occur to no greater extent at the higher temperatures 
than it would at temperatures approaching room temperature. 

Perhaps it should be pointed out again that it is not the hydrogen ion 
activity or pH of the system at elevated temperatures and pressures that is 
quenched and measured, Rather it is the total molar acidity of the system that 
is quenched, and the measurement of pH at room temperature is the means of 
establishing that value. Hydrogen ion and potassium ion activities at the high- 
er temperature conditions of the experiments are estimated from thermody- 
namic considerations. 

Pressure Experiments In addition to the studies at 15000 psi solution 
pressure, experiments were also conducted at higher and lower pressures using 
the 2m KCl solution, The pressure involved is of course essentially water pres- 
sure. The results are given in figure 4, plotted in reference to the 15000 psi 


EFFECT F PRESSURE ON THE EQUILIBRIA 
Mica Pyrophyillite — Andalusite 


Mica — Pyrophyllite — Boehmite 


Mica Kaolinite 


\ 


\ 


\ \ \e 


35000 ps 
5000 
2000 


Displacement of the equilibria with change in pressure, The points in the 
legend are for 2m KCl plotted in reference to the curves at 15000 psi solution pressure. 
The values at the highest temperatures are for 0.65m KCI at 15000 psi solution pressure. 
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curves. The effects of pressure are small, but apparently real. Higher pressures 
increase alteration, lower pressures result in the precipitation of K-feldspar or 
mica for the two reactions respectively, It is evident from the data that com- 


position of the aqueous phase and temperature are more important than solu- 
tion pressure as controls on stability relations. In addition to the experiments 
at high and low pressures, high temperature points for the 15000 psi curves 


are also shown in figure 4, At very high temperatures the 15000 psi curves are 
seen to show a break to higher acidity analogous to that of lower pressure 
curves at lower temperatures. Thus a sheaf of curves corresponding to different 
pressures and fanning out at high temperatures could be illustrated, This rapid 
shift to higher acidity at higher temperatures corresponds to a widening of 
the K-feldspar field relative to the other two. Alteration of K-feldspar would 
occur with diminished intensity with increasing temperature, and at very high 
temperature feldspar would be stable at very high acid concentrations, 
Experiments with Excess Silica.—The role of silica in the displacement of 
the mica-kaolinite equilibrium curve was also investigated, The results are 
shown in figure 5, At 300°C there was no effect on the position of the equilib- 


EFFECT OF EXCESS SILICA 
on the 
Mica—Kaolinite and Mica —Pyroph.—Boehmite 


Equilibria 


L i 
io” io” 


MHCI 


Fig. 5. Experiments with added silica. The points at 400°C are for 2m KCI and the 
value at 300°C is for 4m KCl. These are plotted in reference to the curves at 15000 psi 
solution pressure and silica deficiency. 


251 
500 
& 
400 
300 \ 
200 
= 


252 J. Julian Hemley—-Some Mineralogical Equilibria 


rium, even with the use of highly reactive silicic acid as the source of silica. 
This indicates that mica-kaolinite-quartz is a stable assemblage under these 
conditions. At 400°C, using either cristobalite or silicic acid, there was a small, 
but distinct shift in the equilibrium point. This corresponds to the silication of 
the high alumina phase, boehmite. Here the equilibrium assemblage with ex- 
cess silica is mica-pyrophyllite-cristobalite. Experiments with quartz rather 
than cristobalite were not carried out. The reaction at 400°C would occur with 
somewhat smaller free energy 

At high temperatures the assemblage mica-aluminosilicate-quartz would 
become stable with the breakdown of pyrophyllite. This reaction was not 
demonstrated in the experimentation, although one run at 550°C apparently 
produced some cristobalite. as well as the mica-andalusite-pyrophyllite as- 
semblage, from kaolinite plus pyrophyllite starting materials, However, the 
cristobalite may have been metastably produced in the rapid decomposition of 
kaolinite at this high temperature. 


DISCUSSION OF RESULTS 

Evaluation of the experimental results involves consideration of possible 
errors of experimentation and problems involved in the estimation of equilib- 
rium constants and other thermodynamic parameters of the reactions, Only 
approximate values of these quantities can be obtained in this system because 
of the lack of pertinent data at elevated temperatures and pressures, but some 
thermodynamic conclusions are nevertheless desirable. With such information 
mineral stability relations can be framed in more general terms and the 
several equilibria can be more quantitatively related to other chemical reac- 
tions of geological interest. 

The molar KCl, HCI quotient is an approximate expression of the equilib- 
rium constants of the several reactions studied, For the K-feldspar-mica-silica 
and mica-pyrophyllite-boehmite reactions: 


>KAISi,O, + H+ + 3Si0. + K 


KAI. AISi,0,,(OH) H rAl Si,O,,(OH). 4 {ALO + 


Mice 


K , — , where y+ is the mean activity coefficient of 


Min 
the electrolyte.” Values of the experimental equilibrium quotients are illustrated 
in table 1. They are surprisingly constant in view of experimental difficulties. 
| 
Equilibrium Quotients at 400°C & 15000 psi 
K-feld.-Mica-SiO Mica-Pyroph.-Boeh. 

im KCI 10 

2m KC] 10' 10 

0.65m KCl 10 10 


molal con ration unit indicated is consistent with the further discussion of the 


— 
us | 
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Experimental Errors 

The consistency of the values shown for the experimental equilibrium 
quotients imply that the values are correct to a fairly high degree of accuracy, 
and that the activity coefficient ratios remain essentially constant with mod- 
erately large changes in electrolyte concentration. If the ratios were not con- 
stant and if gross experimental errors were involved in the measurements, for 
example errors as large or larger than the differences in salt concentrations in- 
volved, it would be required that these errors operate in precisely such a way 
as to bring the equilibrium quotients into close agreement. Although this is a 
possibility it is clearly the less plausible view. Such arguments cannot be con- 
sidered proof of a high degree of accuracy of the experimental results, but 
they do constitute a strong indication that the experimental uncertainty is 
small. As observed earlier, an uncertainty of no more than 10°? was indi- 
cated in the reproducibility of points on the equilibrium curves, Errors arising 
from problems in quenching mentioned earlier, which should be the major 
source of uncertainty. would thus appear to be essentially of this magnitude. 

The accuracy of the experimental results can be examined using other 
experimental techniques. For example, a small amount of mineral charge 
might be equilibrated with a large excess of solutions of given K*+/H* ratio, 
and the ratio observed at constant temperature and pressure at which the 
products of the reaction change. Also, for the K-feldspar-mica-silica assemblage 
mica and feldspar can be equilibrated in the manner of the present investiga- 
tion, but without the presence of silica. Fixation of potash does not then take 
place and the acid concentration cannot increase above the initial concentra- 
tion in the charge. Hydrolysis toward higher pH is the only change that may 
occur, Such studies made on the K-feldspar-mica-silica assemblage using 4m 
KCI at 400°C indicated that the maximum quenching error to be expected is 
about 0.5 pH units. The actual error would probably be considerably less, This 
agrees with conclusions already expressed above. 


Activity Coefficients 

Practically no information is available in the chemical literature on the 
activity coeflicients of electrolytes at elevated temperatures and pressures, Cal- 
culation of equilibrium constants must therefore depend upon extrapolation of 
activity coefficients to the higher temperature pressure conditions of the ex- 
periments. 

Temperature Extrapolation.—This point will be discussed, first from the 
standpoint of activity coefficients of the ions based on thermal data at room 
temperature, followed by consideration of the degree of dissociation of the 


electrolytes at elevated temperatures and pressures, The overall activity co- 
efficient will include both of these effects. 
The temperature dependence of the activity coefficient is given by the 


relation: 


Aln y+ H—H 
oT P.M 


H is the partial molal enthalpy of the solute, H® is the partial molal enthalpy 
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at infinite dilution, and L is termed the relative partial molal enthalpy of the 
solute. The mean activity coefficient y+ refers to the ions in molal concentra- 
tion units. The solutions used in the experimentation were made up on a molat 
or volumetric basis, but for the purposes of these approximate calculations, 
and considering that it is concentration ratios that are involved, these can be 
taken as molal values. 

Partial molal quantities for KCl and HCl are known at 25°C and have 
been tabulated by Harned and Owen (1958). Because partial molal heat 
capacities of these electrolytes are unknown at higher temperatures and pres- 
sures extrapolation, of the relative partial molal enthalpy must be based upon 
heat capacity values at low temperature, The relative partial molal enthalpy 
then becomes 


ACpoose (T — 298) 


where ACp is the relative partial molal heat capacity, Incorporating this into 


the above equation, and assuming ACp constant, we obtain the relation: 


Loose — (298) ACp log T 


log 
(2) (2.303) R17 2K 


From knowledge of the activity coefficient and relative enthalpy and heat 
capacity values at some reference temperature such as 298° K the constant | 
can be evaluated. 

Calculations of activity coefficients using this relation are based on the 
assumption that the relative partial molal heat capacity ACp remains constant 
with increase in temperature. In cases where ACp is large, however, this im- 
plies a considerable increase in L and decrease in the activity coefficients at 
higher temperatures. For 2M KCl L 146 cal/mole and ACp 23.8 cal 
mole/deg. at 25°C, At 100°C, L becomes 1790 cal/mole, and continues to in- 
crease to very large values at higher temperatures. The quantity L is a heat of 
dilution with sign changed. Thus if enthalpy changes are to be limited to 
realistic values SCp must decrease with temperature. Lewis and Randall 
(1923, p. 107) have discussed these points and have suggested that at con- 
siderably higher temperatures \Cp 0 is a better basis of approximation 
than the assumption of a constant ACp. The actual trend of the activity co- 
efficients would probably follow some intermediate course between values 


based on Af p const int and Af 1 0. 


Phe properties of mixed electrolytes rather than individual electrolytes 


must be considered in a discussion of the experimental solutions, Not only 
activity coefficients but a relative molal enthalpy and heat capacity values 
of electrolytes lectrolyte mixtures are found to be a function of the total 


concentration or ionic streneth of the solution, The result is that the activity 
coethcient of a minot component in a mixture of ele trolytes tends to follow 
that of the major component 


Mixtures of KCl and HC] at controlled ionic strength have been studied 
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in some detail. An important observation is that at a given concentration of 
KCl, y+-uci is independent of HCl concentration at low acid concentrations. 
Thus the range of experimental pH measurements for a given KCI concentra- 
tion would all involve the same activity coefficient of HCl, The activity coefh- 
cient ratios for the extrapolated activity coefficients are brought much closer 
to unity when the mixed electrolyte, dilute HCl in more concentrated solutions 
of KCl, is considered. This behavior helps explain the approximate constancy 
with change in salt concentration of the equilibrium quotients in the experi- 
mental solutions. 

Extrapolated values from 25°C of the mean activity coefficients of KCl 
and HCl are given in table 2. For the mixed electrolyte the mean activity co- 
efficients for 0.01M HCl in 2M KCI are presented in table 3. These values were 
calculated from experimental data available from 0° to 50°C (Harned and 
Owen, 1958). Such data is also given for KCI solutions, For the mean activity 
coeflicient ratio, the values for 2M KC1/0.01M HCI represent a fair average 
for the experimentation as a whole. These values are presented in table 4. 


TABLE 2 
Mean Activity Coefficients 


ACp constant 


KCl 

0.1M 0.01M 0.001M iM 2M 0.7M 
0.754 0.883 0.982 0.570 0.520 0.579 
0.697 0.850 0.975 0.347 0.372 0.461 
0.646 0.823 0.966 0.204 0.238 0.352 
0.600 0.794 0.960 0.096 0.151 0.261 

202 71 22.6 782 +416 ~85 
2.9 1.4 0.3 pF 23.8 14.0 

ACp 

0.890 0.961 0.661 0.662 0.637 
0.879 0.957 0.740 0.664 0.644 
0.956 0.794 0.691 0.650 
0.956 0.836 0.713 0.653 


al mole 


TABLE 3 
Mean Activity Coefficients of HCl 
0.01M HCl in 2M KCl 
from 25°C from 50°¢ 
ACp const ACp oO AUp const. ACp 0 
0.68 0.71 0.62 0.66 
0.51 0.66 0.32 0.56 
0.40 0.63 0.16 0.51 
0.30 0.61 0.08 0.47 
1120, ACps 1] 


100 

200 

LOO 

ACp 

100 

200 

100 

200 

LOO 
14 34 
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4 
Mean Activity Coeflicient Ratios 


KCI HCI, 0.01M HCl in 2M KCl 


from 50°¢ 
ACp = const. 
100 id $7 0.80 
200 73 1.0 
300 0.60 1.3 
100 0.50 1.4 


Pressure Effect.-\n the above calculations the effect of pressure on the 


activity coefficients of HCl and KCI was neglected. The pressure variation is 


( 


where v° is the partial molal volume of the solute in the standard state, The 


given by: 


conventional standard state would be a hypothetical solution of unit activity 
at | atm. pressure, which obeys the limiting law for infinite dilution, It is ap- 
parent that the effect of pressure on activity coefficients should be relatively 
small. This is fortunate since no data are available at the elevated temperature 
and pressure conditions of the present experiments for partial volumes and 
partial compressibilities of electrolytes. It should be noted that the pressure 
variation of the activity coeflicient is not measured by the difference between 
the partial molal volume of the solute in solution at 1000 atm, versus its partial 
molal volume in the standard state at 1 atm. Rather it is a function of the 
change in the difference between the two volumes at 1000 atm, and 1 atm. 
respectively. As with temperature, changes in activity coefficients of HCl and 
KCl with pressure would parallel one another and tend to cancel, The overall 
effect of pressure on these equilibria is a function of activity coefficients, 
molar volumes and partial molal volumes of the reactant and product materials, 
Although difficult to assess, the uncertainty here is probably considerably less 
than that involved in the temperature extrapolations. 

At lower temperatures, 200° and 300°C, the effect of pressure on activity 
coeflicients would tend to be similar to the effect at room temperature as il- 
lustrated in the table below (from Harned and Owen, p. 507). 


atm 0.1IM HCl IM KCI 
l 0.797 0.6005 


1000 0.802 0.637 


The change at 25°C is negligible, but it is noteworthy that the increase in y+ 
HCl is considerably less than that of y+KCl, Thus, if the same relationship 
held in the mixed electrolytes the effect of pressure on the activity coefficients 
would be to shift the equilibrium quotient to smaller values. 

EMF Data Information on the activity coefficients of pure dilute HCl 


solutions at higher temperatures and moderate pressures has been presented 


) 0 
0.84 
0.95 
1.0 
1.1 
\ \ 
2RT 
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by Lietzke (1955). Lietzke calculated y+ values for 0.017M HCl up to 200°C 
by the methods of partial molal heat content and the Debye-Hiickel theory. 
These two methods agree to within 0.6 percent up to 200°C, Lietzke used these 


values to calculate a theoretical EMF of the hydrogen-silver chloride electrode, 


f’2 H, 


These EMF values were then compared with experimental values of other 
investigators. The experimental potentials were lower than the theoretical 
values by only a few millivolts up to about 150°C. Above 175°C the discrepancy 
increased, and at 200°C and 42 atm. the theoretical potential was 0.489 V, the 
experimental 0.469 V. Above this temperature both potentials fall off rapidly 
because of rapid decrease in the partial fugacity of H.. Higher pressure of 
course increases the potential at any temperature, the maximum pressure used 
being approximately 42 atm. 

The discrepancy at 200°C would correspond to a difference in the ob- 
served and theoretical activity coefficients of only about 5 percent. It is in the 
wrong direction to be attributed to a decreasing activity coefficient of HCl. 
The work of Lietzke is therefore evidence that HCl remains a strong electrolyte 
at these moderately high temperature pressure conditions. This would support 
the assumption that the activity coefficient values listed in tables 2 and 3 are 
approximately correct. In the mixed electrolytes HC] would have the greater 
uncertainty because of the rapidity with which the activity coefficients of HCl 
are decreasing with temperature at low temperatures in these solutions, 

Solubility Data—Some further insight into the activity coefficients of 
KCI solutions is afforded by the work of Keevil (1942) on vapor pressures 
and solubilities in the KCI-H.O system at elevated temperatures and pressures. 
The halides of sodium and potassium fall into the class of compounds whose 
solubility curves are continuous to the melting point of the salt. Very great 
aqueous solubilities are therefore realized, provided the external pressure is 
maintained greater than the vapor pressure. 

The activity coefficient for the saturated solution can be calculated, The 
free energy of fusion at any temperature and the free energy of solution can 
be equated since their sum must be zero. 


KCl.) a.) AF®° (fusion) 
KCl aig.) KClisom) AF = RT Ina 
KCl.) AF 
AF AHo — ACp Tin + IT RT In a 
and AH (fusion) AHo + ACpT. 
lhe constants AHo and I can be evaluated at the melting point (1043°K) 
from the heat of fusion, and the difference in heat capacity of solid and liquid 


KCl. Activities may then be computed, and the activity coefficients correspond- 
ing to the solubility data of Keevil obtained. 


+H, + AgCl = H+ + Cl- + Ag 
RT (M y+)? 
E E° In 
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For better comparison with the extrapolated activity coefficient values 
shown in tables 2 and 3, activity coeflicients derived from Keevil’s data may 
be computed on the basis of a moles 1000g. of solution standard state rather 
than the pure liquid standard state. The total free energy change is then given 


by Al + RT In m? y+? 0. where the second 
m” Ke 


term involves the change in standard state and the third term refers to Keevil’s 
solutions. Using the thermal data of Kelley (1949), the mean activity coefh- 
cient of KCI at 400°C is then 0.72. At 500°C it is 0.81, Although the standard 
state now used is not the same as that of the extrapolated activity coeflicients, 
a hypothetical 1 molal solution obeying the limiting law for infinite dilution, 
these two standard state fugacities and corresponding activity coefficients will 
not differ greatly from one another if the solutions are not far from ideal. Also. 
both will be close to unity. The value of 0.72 at 400°C tends to indicate such 
behavior 

Finally it will be noted that the activity coefficients given by Keevil’s data 
refer to the saturated solution at the saturation vapor pressure (e.g, 125 atm. 
at 400°C), At 15000 psi the solubility should be somewhat higher and the 
activity coeflicient somewhat lower. Also, it is possible that at much lower 
concentrations activity coeflicients would be somewhat smaller since y+ values 
typically pass through a minimum at some moderately low concentration and 
then rise to higher values. 

These computations constitute an approximate check on the accuracy of 
the extrapolated activity coeflicients given in table 2, It is apparent from these 
results that the values in table 2 imply a mean activity coeflicient of KCI that 
is certainly of the correct order of magnitude. Furthermore, although the ab- 
solute values of the activity coefficients of KC] and HCl are subject to con- 
siderable uncertainty, the uncertainty in the ratio of the activity coefficients 
would be muc h smaller, The values should compare fairly closely to those in- 
dicated in table 4. due to the behavior of the mixed electrolytes. 

Conductance Data.—The strength of HCl and KCl as electrolytes at high 
temperatures and pressures has been evaluated by Frank (1956). Conduc- 
tances, A, were determined at temperatures from 250° to 750°C and at pres- 
sures up to 2500 atm. Using the method of Shedlovsky. Frank evaluated A,,. 
the limiting equivalent conductance in dilute solution. The degree of dissocia- 


tion of the electrolyte is represented by the relation a 5 . The equilib- 


rium dissociation constant K . where c is the concentration, The 


l—a 


activity coefficients for the dilute solutions involved in such studies can be 
estimated from the Debye-Hiickel equation. Dissociation constants obtained by 
Frank are given in table 5 

The most important control on the strength of the electrolytes is the steam 
density. At densities greater than 0.8 2/cm* at 400°C. for example, « for KCl, 
HCl and KOH eradually approac h unity, but at lower densities @ decreases 
rapidly, HCl being affected to the greatest extent. At constant steam density 
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a decreases slowly with increasing temperature. At low temperatures densities 
are high and @ is essentially unity. These findings indicate that HCl, KCl and 
KOH remain quite strong electrolytes at elevated temperatures and pressures 


provided the water density remains high. These general results are consistent 
with the work of Lietzke in the lower range of temperatures. 

Possible values for the mean activity coefficient of 2M KCl and 0.01M 
HCL at 400°C would be 0.3 and 0.1 respectively. These are taken from the 
data presented in tables 2 and 3, which, as already noted, are substantiated by 
the experimental studies that have been discussed, The effect of non-dissocia- 
tion on these extrapolated values can be computed from the ionization con- 
stants in table 5. Assume the activity coeflicients of the undissociated species 
to be 1. 


aCl 
mK Cl ak 
0.3 


Ke 


aCl 0.20 
aH* aCl 
mKCl aH 
0.1] 


a~KCl ak* aCl 
a~HCl aCl 


The total mean activity coeflicient ratio is , °.°° whereas the ex- 
trapolated ratio is 


TABLE 5 
Dissociation Constants 
KCl 
density g/cm 


0.3 0.5 0.6 
21-10 
56-10 2.75-10°° 
6.0-10°° 
10.3-10 3: 15-10 
2.9-10° 
6.3-10 . 1.00-10 


HCl 
100 1.0-10 11-10° 7 7.0-:10°° 
00 7-10°° t 6.3-10° 29.10 
600 1.8-10 29.10 1.0-10° 
700 0.7-10 3 1.4-10 
KOH 
100 90-10" 2. 2.7-10 5-107? 15-107 
500 1.9-10°° 1.9-10 107 3.4-10 
600 0.8-10 We 1.3-10°° 
700 0.4-10 2 1.0-10 


after Frank 


ak 
7.0-10 
aH 2.58 LO 
102.59 
107-89 
( 0.7 0.8. 
100 31-1077 
150 
12-10% 
950 
600 
650 
750 1.2-10-* 
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The above computations serve to illustrate the extent to which non- 
dissociation of the electrolytes might contribute to the overall activity coefh- 
cient ratio. The effect would be significant at 400°C, but would become 
relatively less important at lower temperatures. 

Conclusions._In summary, on the basis of the preceding discussions it 
appears that the experimental equilibrium quotients would deviate from true 
equilibrium constants of the reactions by no more than a factor of 5 up to 
temperatures of 350°C and possibly higher. Activity coefficient ratios given in 
table 4 suggest an uncertainty no greater than a factor of , A but these values 
may not be reliable to that extent. Because of the relatively rapid change in 
the mean activity coefhicients of HCl from 25° to 50°C there is considerable 
uncertainty involved in the extrapolation. Electrode potential measurements to 
only LOO° to 150°C would probably do much to eliminate this uncertainty, 


Pressure rperiments 
In figure 4 are indicated the results of experiments in which the pressure 
was varied from 15000 psi for the 2M KCI solution. The points are plotted in 
reference to the 15000 psi curves, The effect on the equilibrium quotient is 
small 


Phe free energy change occasioned by a change in pressure is given by 


the relation AVdp. where AV is the net volume change of the reactants 


and produc ry From this the shift in the equilibrium quotient can be deduced. 
The partial molal volumes of the mixed elec trolytes are probably very similar. 
For the individual electrolytes the difference in partial volumes of 0.01M HCl 
and 2ML KCI at 25°C is about 13 cc. and in the mixture this difference should 
become even smaller. Also, the density of the experimental solutions at 400°C 
and 1000 atm. would be quite large. slightly greater than 0.7 ¢/cm*, and large 
changes in partial molal volumes would not be expected, For the solids the 
volume change can be estimated from the known densities of the natural min- 
erals, bearing in mind that the densities of the imperfectly crystalline synthetic 
material is probably somewhat less than ideal. For the K-feldspar-mica-silica 
hydrolysis reaction the molar volume change is about —20 cc/mole H*, The 
shift in the equilibrium quotient would be about 10°-'* for a pressure increase 
of 1000 atm, This is approximately the observed change, This suggests that 
differences in partial molal volumes and changes in the a¢ tivity coeflicient ratio 
for the mixed electrolytes are small, although some compensating changes in 
activity coefficients and molar and partial molal volumes might also occur. 

\ change to lowe pressure produces a somewhat greater observed effect. 
This is consistent with the pattern of the dissociation constants noted in table 5. 

Phe effect of pressure on activity coefficients of HCl and KCl at low tem- 
peratures was noted on a previous page. Although the effect is small, the 
activity coeflicient of KCl is increased considerably more than that of HCl. 
Similarly, it will be observed that the high pressure points on the kaolinite- 
mica curve at 300°C tend to shift the curve slightly to the left. If this effect is 
real it is a further indication that extrapolation of activity coefficient ratios to 
at least 300°C involves little error. 
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A marked displacement of the 15000 psi curve at high temperatures is 
also shown in figure 4, Although this shift to higher acidity illustrates the 
same temperature effect as that evident in the lower pressure curves at lower 
temperature, the process is probably somewhat -different. The rapid rate of 
change probably results from passing into a two phase critical region, Parti- 
tion of more HCl in the lighter phase, with correspondingly lower activity co- 
eflicient, produces the observed effect. 


Experiments with Excess Silica 

The influence of added silica on the mica-kaolinite and mica-pyrophyllite- 
boehmite equilibria is illustrated in figure 5. The point at 300°C shows no 
effect on the equilibrium quotient. Silica apparently takes no part in the net 
reaction, which indicates that kaolinite plus quartz would be a stable assem- 
blage at this temperature. The absence of reaction between kaolinite and 
quartz at 300°C, and at temperatures in excess of 350°C, in experiments of 4 
months duration at 15000 psi water pressure is additional evidence of such a 
stability relation. 

At 400°C boehmite is silicated to pyrophyllite, which shifts the equilib- 
rium quotient to a new value corresponding to the assemblage mica-pyrophy]l- 
lite-cristobalite. The displacement is a measure of the free energy of the 
reaction 


Al.O.-H.O + 4Si0.;c,) 


The AF 925 cal, mole, with a large uncertainty because of the small dis- 
placement. If the experiments had been continued at higher temperature there 
would presumably have been similar silication of andalusite to pyrophyllite. 

At higher temperatures the assemblage mica-aluminosilicate-quartz would 
become stable with the breakdown of pyrophyllite. The decomposition of 
pyrophyllite occurs at 575°C at 15000 psi water pressure (Roy and Osborn, 
1954) (Kennedy, 1955). The product is mullite or mullite-sillimanite solid 
solution, At still higher temperatures, the assemblage K-feldspar-aluminosili- 
cate-quartz becomes stable with the thermal decomposition of mica, as shown 
by Yoder and Eugster (1955). 


Activity of Water 
The foregoing discussion has been concerned only with those equilibria 


in which water takes no part in the net reaction, But in the decomposition of 
mica to kaolinite: 


K a=Kec1 Mxe 


Although the experimental solutions were moderately concentrated with 
respect to KCl in the sense of molar concentration, the mole fraction of water 
was nevertheless greatly preponderant over that of KCl. (e.g. Nxo,) = 0.07 in 
4m KCl). This of course means that deviation from ideality on the part of the 
solvent water would be relatively much smaller than that of the solute, and it 
is for this reason that activities and activity coefficients of solvents, even in 
fairly concentrated solutions, are not far from unity. 
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An approximate measure of the departure of au.o from unity is obtained 
in terms of general thermodynamic considerations. The standard free energy 
of the mica-kaolinite reaction corresponds to the free energy change when 
water is in its standard state, at unity activity. The experimental solutions at 
15000 psi consisted of a single homogeneous phase of liquid character, and 
therefore the liquid state would be the logical choice for the standard state, By 
convention the standard state would be pure liquid water at 1 atmosphere 
pressure and the respective temperature. In other words, the standard state 
would be a hypothetical state of pure, super-heated liquid water at 1 atmos- 
phere pressure. The change in activity of water in the solutions can be viewed 
as a combination of two processes: The change in free energy when pure water 
goes from 1 atm. to 1000 atm. plus the change in free energy at 1000 atm. 
when pure water goes to watet in the elec trolyte solution. For the first process 


AF is given by { vdp, where v is a function of p 


AF = RT Ina RT In vdp 
e 
The change in molar volume v would not be great for liquid water in the 
two states. This can be seen in the P-V-T data of compressed water (Dorsey, 
1940). Assuming an average specific volume of 2, which is certainly too large 


f 
a value, for the pressure interval 1 to 1000 atm, at 300°C an,o 22. 


f 
Thus by virtue of the increase in pressure, the fugacity would be increased by 
a factor of no more than 2 over the standard state fugacity. In the second 


proce ess: 


(soln) 


AF = RT In = RT In 


H2° (pure) 


The sum of these free energies would then correspond to the activity of water 
in the solution at 1000 atm. referred to the proper standard state. 
soln fi (soln) 
Al RT Ina ‘ ) RT In f 
1140 

The decrease in fugacity or activity occasioned by the transport of water from 
the pure state to the solution would be small. From these considerations it 
would be expected that the activity of water in the experimental solutions 
would be less than 2, and possibly slightly less than 1 in the more concentrated 
solutions. 

The above discussion could have been based upon any standard state for 
water, Plausible alternatives would have been pure gaseous water at 1 atm. 
pressure or water at 1000 atm. pressure. If the gas standard state were chosen, 
the activity change from H.O (2. 1 atm.) to H.O (2, 1000 atm.) could be 
calculated from known fugacity coeflicients of water (Holser, 1954). 

Consideration of au.o helps explain the difference in curvature of the 


equilibrium constant curves below 350°C, Figure 2 is a plot of the K values 
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in which the au.o for the mica-kaolinite reaction is assumed to be 1. The latter 
curve appears to be essentially the metastable extension of the mica-pyrophyl- 
lite-boehmite curve below 350°C, and suggests that the true curve should have 
a break in slope and trend somewhat to the right, Holser’s data indicate that 
the fugacity of water at 200°C and 1000 atm. is only 21, as compared to 109 
at 300°C. Thus, when relatively smaller activities of water are substituted into 
the expression for the equilibrium constant at lower temperatures, the constant 
is increased and the curve shifted relatively to the right. At the transition tem- 
perature, unless the activity of water happened to be unity, the two curves 
would be somewhat displaced from one another. The dashed curve illustrated 
would correspond to anyzo 0.87 at 200°C and axn.0 = 1 at 350°C, Exact 
parallelism of the curves would not be expected, since that would imply identi- 
cal heats of reaction. 


Free Energy and Heat of Reaction 
Values for the standard heats of reaction for the results at 15000 psi are 
given by the curves shown in figure 6, These are derived from the relation 


dln K AH° 
oT |] P.M RT? 


(300°) 


(350°) 


= iL 
10! 10? 10 


Equilibrium Quotient 


Fig. 6. Plot of the equilibrium quotient versus the reciprocal of the absolute tem- 
perature. 
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The heats of reaction are exothermic, as would be expected for hydration or 
hydrolysis reactions, The units are kcal/mole H+, These results are based on 
the assumption that the equilibrium quotient is a close approximation of the 
equilibrium constant of the reactions as already discussed. The uncertainty in- 


dicated would correspond to an experimental reproducibility of 10*°* in the 


equilibrium constant. 

Although the mica-kaolinite curve is plotted on the assumption of au,o 
1, a small break in slope is indicated at the kaolinite decomposition tempera- 
ture. A more pronounced change in slope, which for example would cor- 
respond to the dashed curve of figure 2, is probably correct. This would imply 
a more negative heat of reaction. 

Standard free energies of reaction are tabulated in table 6 as given by 
the relation 

AF RT Ink 

The free energy drive for hydrolysis decreases with increase in temperature. 
The standard entropy change is obtained from the free energy and enthalpy 
values. Average values of —11.3 and —14.8 cal/mole/deg. are indicated for the 
hydrolysis of K-feldspat and mica respectively. The dashed mica-kaolinite 
curve would cor respond to a more negative entropy of reaction. 


TABLE 6 
Standard Free Energy of Reaction 


keal/mole H 


300 100 
K-feld.-Mica-SiO, 0.43 9.30 + 0.52 3.31 + 0.61 
Mica-Kaolin; et 70 = 0.43 5.24 + 0.52 1.00 + 0.61 


At 25°C the standard entropy change for the hydrolysis reactions can be 
estimated by procedures such as those described by Latimer (1952). The 
values obtained, however, are rather uncertain in the case of the present re- 
actions involving hydrated silicate minerals. A point that might be mentioned 
is that the standard entropy change for these reactions at room temperature is 
positive. This positive change in entropy becomes negative at higher tempera- 
tures, This is accomplished by the relative increase in ACp for HCI in the 
mixed electrolyte as the temperature increases. 

Thermodynamic relations between K-feldspar, mica and kaolinite at 25°C 
have been considered by Garrels (1957). On the basis of field relations and 
hydrolysis experiments on K-feldspar and mica, equilibrium constants of 107° 
and 10°° were proposed for the hydrolysis of K-feldspar and K-mica respec- 
tively. These are minimum values, with an estimated uncertainty of about 
10° in the equilibrium constant or 1 kilocalorie in the standard free energy 
of reaction per mole of H*. It can be seen in figure 2, or in an extrapolation 
of the curves for the heat of reaction,’ that the mica-kaolinite curve would ex- 
trapolate to approximately the value proposed by Garrels. The K-feldspar 
hydrolysis curve would extend to a value somewhat higher than that given by 
Garrels. The diser pancy in this instance suggests an error of about 2 kilo- 


500 
7.42 + 0.71 
ee 2.83 + 0.71 
4 
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calories in Garrels’ values. The agreement between Garrels’ values and the 
results of the present experimentation is satisfactory. 

Measurements of activity coefficients at moderately low temperature- 
pressure conditions are very much needed. They would be of great value in- 
asmuch as the rate of change in the thermodynamic properties of these electro- 
lyte solutions is rapid in this low temperature range. 


Phase Relations 

Vineral Phases.—These relations in the system K.O—AI,O,—SiO.—H.O 
are simply stated. In the hydrolysis of K-feldspar to mica plus silica the 
minimum number of components is 4. These may be taken as K-feldspar, mica, 
silica, water or K,O, Al.O,, SiO., and H.O. Applying the phase rule, the 
variance w 14.42 2. At a given temperature, the system will be at 
equilibrium over a range of pressures, or at a given pressure over a range of 
temperatures. If the temperature and pressure are fixed the system will attain 


equilibrium at a unique K*,/H* ratio and unique concentration of KOH. If 


chloride is added to the system a new component is introduced, This will in- 
crease the variance to 3. Here, in order to define the composition of the aqueous 
phase the concentration of either KCl, HCl. or total chloride must be specified, 
even though the K* H®* ratio will be the same as in the previous case at a 
given temperature and pressure, 

The same type of analysis applies to all of the assemblages encountered 
in this study. For mica-kaolinite-water only 3 phases exist, which would imply 
a variance of 3. Here, however, the number of independently variable oxide 
components can be reduced to 3 because of the restriction of the same ratio of 
Si Al in the mica and kaolinite. Addition of acid increases the number of 
components in the same manner as discussed above. In the thermal decomposi- 
tion of muscovite to K-feldspar plus corundum, on the other hand, the system 
has a variance of only 1. A univariant dehydration curve describes this type 
of reaction, but the curve is not unique if added electrolytes are present, 

The activities of Al,O, and SiO.. as well as those of the other components, 
are fixed at a given pressure, temperature and total chloride concentration. 
Addition of alumina to the K-feldspar-mica-silica assemblage, for example, 
would only produce more mica at the constant K+ /H* ratio. If added to the 
mica-kaolinite reaction it would simply precipitate as excess alumina. In this 
connection it is of interest that in the hydrolysis of mica to kaolinite in the 
more saline solutions a small amount of boehmite appears with kaolinite in the 
reaction product, This superimposed reaction results from greater solubility of 
silica under these conditions, but it does not affect the equilibrium constant 
for the reaction under consideration. 

\ point of special interest is the question of equilibrium between as- 
semblages at the thermal breakdown of kaolinite to pyrophyllite plus boehmite. 
Here 5 phases would be in equilibrium and 4-5+2 1. At a given pressure 
there would be only one temperature, the thermodynamic decomposition tem- 
perature of kaolinite, at which all of the phases could exist in equilibrium, In 


* These do not extrapolate as straight lines, but must steepen in slope at lower tempera- 


tures consistent with the entropy change already noted. 
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this instance the activity of water is particularly important, for it determines 
the decomposition temperature of kaolinite at any given salt concentration, In- 
creasing the concentration of dissolved electrolyte decreases the activity of 
water and lowers the decomposition temperature. 

The temperature of decomposition of kaolinite at 15000 psi pure water 
pressure was not established in the present study. Conclusions as to the de- 
composition temperature of kaolinite in the electrolyte solutions are based on 
the points at 352°+5°C shown in figure 1. These results incidentally, are 
consistent with the experimental results of Gruner (1944), The decomposition 
temperatures of kaolinite in the 2m and 4m KCI solutions are indicated only 
diagramatically and may be somewhat lower, especially in the 4m KCI solu- 
tion, Roy and Osborn (1954) and Sand (1956) found that kaolinite decom- 
poses at 400-405°C at 15000 psi water pressure. The product, a material which 
does not occur in nature, was named hydralsite by Roy and Osborn, Pyrophyl- 
lite also was generally present. Hydralsite was not observed in the present ex- 
perimentation using electrolyte solutions. 

KCI—H,,O System.—Phase relations in the KC]—H.O system are im- 
portant to the experimentation. These are described in other papers. Keevil 
(1942) determined the solubility curve and saturation vapor pressures in the 
system KCI—H.O up to 600°C. A paper by Copeland. Silverman and Benson 
(1953) on the system NaCl——H.O describes phase relations at low concentra- 


tions in the supercritical region 


Effect of Composition 

This experimental work has been concerned exclusively with the evalua- 
tion of certain mineralogical equilibria in which pure solid phases in the 
system K.OQ—ALO,—SiO.—H.O are involved. If. as would occur in nature. 
other components are present that would enter the reaction, equilibria would 
he correspondingly changed. For instance, a small amount of Na in the K- 
feldspar would result in a correspondingly small change in the equilibrium 
constant for hydrolysis to mica, and any change in the composition of the 
resulting mica would have a similar effect. The magnitude of an equilibrium 
constant expresses the free ene rey change of a reaction, and therefore any 
change in the free energies of formation of any of the reactant and product 
phases must affect its value. For the hydrolysis of K-feldspar to mica + silica 


Amica AK 


. With only small 


the equilibrium expression now becomes 
QK-teldspar Gu 


changes in the composition of the solids the equilibrium constant would be but 
little affected. 


GEOLOGICAL APPLICATIONS 
In the applications of the experimental results to geological problems the 
ratio of the activity of aqueous potassium to the activity of acid in the hy- 
drolysis reactions for K-feldspar and mica will be referred to simply as the 
K* H®* ratio. From the standpoint of the effect of solution composition on 
stability relations this quantity is the important thermodynamic consideration, 


regardless of what species of acids or potassium compounds may actually be 
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present in the natural system. The equilibrium K*+/H* values are, of course, 
independent of the H* activity of the solution, its acid or basic character, The 
geologic applications of the present experimentation have been, for the most 
part, stated or implied in discussions on previous pages of the several equilib- 
ria and equilibrium curves. Only a review of some of the salient points is ap- 
propriate here. 


Hydrothermal Processes 


In the experimental curves of figure 2 it is observed that the equilibrium 
K* H®* ratios trend toward higher acidities at higher temperatures. This slope 
means that the potential of hydrothermal solution to alter its wall rock as it 


ascends from environments of higher temperature would not necessarily 
diminish as its temperature falls. In fact, that ability would increase if the rate 
of migration were sufficiently rapid in comparison with the rate of reaction 
with the wall rock. In order to avoid confusion one should distinguish two 
different vantage points of analysis (Meyer, 1950), that which transpires at 
any given point in the vein or system with the passage of time, and that which 
takes place in an increment of solution as it rises from depth to higher levels. 
Considering first the solution, three competing processes would operate, (1) 
as the temperature falls the reactive potential of the solutions would increase 
due to the negative slope of the equilibrium curves and, at higher tempera- 
tures, the increase in dissociation of the acid components, (2) increased re- 
action with the wall rocks would tend to compensate for the increased reac- 
tivity of the solutions, (3) lower temperature would result in lower rates of 
reaction. The net effect would be a complex combination of these processes, It 
is at once apparent that the total amount of alteration in a hydrothermal 
system would be small at points representing the two extremes of very high 
and very low temperature, with maximum alteration somewhere in between. 
At any intermediate point conditions would change with time, In detail this 
change would be a function of many factors such as rates of reaction, fluid 
flow. heat transfer and diffusion. even assuming for the moment constant com- 
position of solutions at “source”. The over-all effect, however, would be in- 
creased reaction with the wall rock with falling temperature which would tend 
to maintain the K*,H* ratio at an equilibrum value, Continuous leaching 
and steadily increasing chemical insulation in the walls along the preceding 
path, on the other hand, would tend to produc e a gradual decrease in the ratio. 
But a smaller ratio again means greater reactivity of the solutions, In this way 
the system would constantly seek to establish equilibrium at every point along 
its path, in accordance in fact with the principle of Le Chatelier. 

The above mentioned process of progressive chemical insulation might be 
accompanied by a similar thermal effect. A rising temperature at any level 
would tend to counteract any increased reactive potential of the solutions. 
Similarly, from the pressure experiments it is apparent that decreasing pressure 
would tend to compensate for the effects of decreasing temperature, In figure 
t the displacements of the equilibrium curves by pressure show that an in- 
crease in pressure increases the extent of hydrolysis whereas a decrease in 
pressure shifts the equilibria in the opposite direction, If pressure effects were 


208 J. Julian Hemley—Some Mineralogical Equilibria 


dominant at any point, as would very likely be the case for surface or near- 
surface processes, the result would be deposition rather than alteration, This. 


then, is a logical source of the adularia observed in some hot springs and 
“shallow” precious metal vein deposits. 

The experimental results have considerable application in the interpreta- 
tion of alteration zoning and the development of various alteration patterns. 
These problems. however. have been more fully discussed in another paper 


(Meyer and Hemley. 1957) 


Vetamorphism 

Some implications of the experimental results to rock metamorphism 
might also be noted Metamorphi and hydrothermal effects are, of course. 
gradational. The same chemical stability relations apply to mineral genesis in 
either instance 

The appearance of the assemblage muscovite-chlorite or muscovite-biotite. 
often together with albite epidote. is one of the first changes to occur in the 
development of low grade regional metamorphism. If this transition is con- 
sidered the lower limit of the Greenschist facies. it is clear that its position 
would be a function of the K H> ratio of the environmental solutions as well 
as the pressure-temperature conditions. At higher K*,/H®* ratios the facies 
boundary would shift to lower temperatures. 

At higher temperatures these equilibrium ratios or constants are likewise 
an intrinsic part ol the mineralogical equilibria. In the A mphibolite facies. for 
example. two important general groups of assemblages containing excess silica 
may be designated. These are defined as containing an excess and a deficiency 
respectively of K.O (Turner and Verhoogen, 1951, p. 448). In the former, k- 
fe ldsp ir co-exists with muscovite. whereas in the latter andalusite or andalusite 
plus muscovite is observed. This relationship is consistent with the equilibrium 
curves of figures | and 2. The muscovite-andalusite and muscovite-K-feldspar 
curves are separated, wi the feldspar incompatible with andalusite. The 
metamorphic processes responsible for the development of these two meta- 
morphic assemblages. were therefore characterized by different K* /H* ratios. 
Thus, in addition to the difference in K.O content of the rocks as they now 
appear, the andalusite bearing assemblage could be described as the result of 
processes operating in a lower 1 inge of kK H values, The converse would 
he true of the other assemblage. In this way the investigation of chemical 
equilibria of the type considered in the present study would be helpful in the 


analysis and interpretation of various metamorphic problems. 


FUTURE WORK 

There are a great many possibilties for continued research along these 
lines. These would include fundamental studies of electrolyte solutions at ele- 
vated temperatures and pressures, expansion of the experimental system to in- 
clude more components, and kinetic studies on alteration reactions, 

Mineralogical equilibria with Na,Q—CaO—AI.O,—SiO.—H.O and 
K.O—Meg0—AI.0,—SiO.—H.0O as components of the alteration system are 
currently being studied. The alteration and stability of plagioclase of varied 
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composition as a function of pressure, temperatures and chemical environment 


is of course a problem of major geological importance. Similarly, equilibria 


between possible alteration products such as montmorillonite, paragonite, and 
zeolites are of interest. In the alteration system containing Mg0—K.0, stability 
relations involving phases such as phlogopite, chlorite, and serpentine are 
analogous problems that need to be investigated. 
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THE MAGNETIC FIELD 
IN A PRIMARY METEORITE BODY 
JOHN F. LOVERING 
Dept. of Geophysics, Australian National University, Canberra, Australia 


ABSTRACT. The magnetic properties and mineralogical composition of the Moore 
County eucrite indicate that, while situated within about 10 km of the surface of the 
silicate mantle of a primary meteorite body, this meteorite cooled through the Curie tem- 
perature of its magnetic constituents (about 560°C) while in a magnetic field. If this field 
associated with the primary body was generated within a fluid metallic core, then at this 
time temperatures within the core must have been greater than about 1700°C, Further- 
more the inclination of the direction of magnetization with the horizontal plane of crystal 
layering observed in the meteorite suggests that the Moore County eucrite formed close to 
either the 10°N or 10°S magnetic latitude of the primary body. 
INTRODUCTION 

The most promising explanations of the Earth’s magnetic field invoke the 
action of hydrodynamical motions within its conducting fluid core (Hide, 
1956). As it appears that iron meteorites formed by the solidification of an 
originally fluid metallic core within a primary meteorite body (Lovering, 
1957a) it is not unlikely that a magnetic field was also associated with the 
primary body. 

\ preliminary study was made of the magnetic properties of various pri- 
mary meteorites to see if any evidence could be found of the existence of a 
magnetic field in the primary body. Those primary meteorites containing a 
considerable proportion of iron-nickel metallic phases (i.e. iron meteorites and 
pallasites) were all found to have relatively low remanent magnetic intensities 
(table 1) comparable to those acquired by nickel-iron alloys in the Earth’s 
field. These meteorites have also been shown to be magnetically “soft” (F. D. 
Stacey, personal communication) so that their observed remanent magnetic in- 
tensities were probably acquired in the terrestrial field after they had fallen on 
the Earth’s surface. 

Certain achondritic stony meteorites (i.e, the eucrites or basaltic achon- 
drites) have apparently crystallized in the outer regions of the silicate mantle 
of the primary body. These meteorites have compositions and textures very 
similar to those of certain terrestrial basic igneous rocks. such as basalts and 
particularly dolerites, whose magnetic properties are widely used to determine 
terrestrial palaeomagnetic fields. Probably the most intensively studied of these 
achondrites is the Moore County eucrite (Hess and Henderson. 1949) and it 
was this meteorite which was selected for the present study. 


REMANENT MAGNETIZATION OF THE MOORE COUNTY EUCRITE 

A fragment of the Moore County eucrite was found to have a remanent 
magnetization with an intensity (table 1) similar to that of the lowest values 
observed for Cenozoic basic igneous rocks which have cooled through their 
Curie temperatures in the Earth’s magnetic field. The magnetic constituents of 
the Moore County eucrite make up less than 2.5 percent of the meteorite and 
a preliminary thermomagnetic study indicates that there are two main com- 
ponents contributing to its remanent magnetization. The magnetically most im- 
portant phase is magnetite with a Curie temperature of about 560°C while 


there are minor metallic phases which are probably demagnetized or only 
weakly magnetic by analogy with the iron meteorites. 


The most likely way in which the Moore County eucrite could have ac- 
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Taste 
Remanent Magnetic Intensities of Some Primary Meteorites 
Remanent 
Magnetic 
Nickel content Intensity 


of Metal Phase*  (10°°c.g.s. 
Py pe Name Classification (%) per ec.) 


Mount Stirlir Coarse octahedrite 6.93 16.1 


Iron Kyancutta Medium octahedrite 8.28 29.0 
Meteorites Mungindi No. 2 Fine octahedrite 12.36 60.3 
Catherina Nickel-rich ataxite 33.60* $39.0 


Stony-lron Pallasite (metal 


9.20 
Meteorites phase only) sg 


Stony 
(Achondritic) Moore County 
Meteorites 


Lovering, et al. (1957) This work. 


quired its observed remanent magnetization is by cooling through the Curie 
temperature of its constituent magnetite (560°C) while situated in a magnetic 
field and presumably while still forming part of the outer mantle of the pri- 
mary body. An important corollary is that while the metal core of the primary 
body was at some temperature above that of the lowest melting point iron- 
nickel alloy within the core |i.e. greater than about 1700°C, assuming maxi- 
mum core pressures of about 10* to 10° atmospheres (Lovering, 1957b) | and 
thus still at least partially fluid and capable of field producing motions, the 
outer silicate mantle had crystallized and cooled to a temperature of below 
960°C. This conclusion provides important data concerning the cooling history 
of the primary body as well as its internal temperature distribution. However 
it does not agree with the conclusions of a study of the Moore County eucrite 
by Hess and Henderson (1949), They found that the lowest temperature 
reached by the meteorite in the mantle of the primary body and immediately 
prior to its disruption was about 1135 + 65°C, This figure was based on two 
observations; firstly, the limited (about 11 percent) inversion of the pigeonitic 


pyroxene to hypersthene and secondly the occurrence of tridymite in the 


Moore County. However it is well established that the presence of impurity ions 


can buttress the tridymite lattice and prevent its transformation to quartz when 
the temperature is lowered (Eitel, 1954). Since the tridymite in the Moore 
County meteorite is interstitial and evidently the last of the essential constitu- 


ents to crystallize, it may well contain a relatively high content of impurity 
ions. In point of fact Hess and Henderson (1949) have reported that the tridy- 
mite in the Moore County meteorite is crowded with minute inclusions, indicat- 
ing a high impurity content in the silica. Thus the presence of tridymite in the 


\Moore County eucrite gives no unequivocal data concerning the cooling history 


of the meteorite. However it does tell us something of the pressures operating 
on this meteorite in th mantle ol the primary hody. It has been calculated 
from thermodynamic dat \losesman and Pitzer, 1941)! that the stability field 


diagram for SiOz has recently been pub 
(1958) which is considerably different from the 
(1941) on the basis of thermodynamic calcula 
that tridvmite is stable up to pressures of about 
heres given by the theoretical calculations, Tuttle 
s not possible to decide which diagram is correct 

sults it would appear that the Moore County 
1) km in the mantle of a primary meteorite 
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of tridymite disappears at pressures above 650 atmospheres. Above this pressure 
liquid silica crystallizes to cristobalite and above 1, 160 atmospheres cristobalite 
is no longer stable and only quartz crystallizes. If the primary meteorite body 
had interior pressures of between 10‘ and 10° atmospheres (Lovering, 1957b) 
then a pressure of 650 atmospheres would correspond to a depth of about 10 
km in the mantle of such a body, The existence of tridymite in the Moore 
County eucrite suggests that this meteorite formed within about 10 km of the 
surface of the primary body, provided that impurities have little effect on the 
pressure above which tridymite is unstable. 

The evidence of the apparently restricted inversion of the pigeonite in the 
Moore County eucrite is still considered by Dr. H. H. Hess (personal com- 
munication) to be indicative of rapid chilling from a temperature of about 
1135°C after it became a meteorite. Brown (1957) has studied the pyroxenes 
of the Skaergaard intrusion and has shown how uninverted and partially in- 
verted pigeonite need not have been due to rapid cooling of the magma, but 
rather to the sluggishness of the inversion of the intermediate and ferri-ferrous 
pigeonites at a stage when an appreciable time interval separates the tempera- 
tures of crystallization and inversion. The Moore County eucrite has already 
been shown to have formed in the outer regions of the mantle of the primary 
body and its partially inverted pigeonite has a composition Ca; Mgy;.5 Fes;.s 
(Hess and Henderson, 1949). From Brown's conclusions it would not seem at 
all unlikely that the pigeonite of the Moore County eucrite could have re- 
mained only partially inverted even after cooling to temperatures below 500°C 
before the primary body was disrupted. 


SITE OF CRYSTALLIZATION OF THE MOORE COUNTY EUCRITE 


Hess and Henderson (1949) have reported that the pyroxene and plagio- 
clase crystals in the Moore County eucrite show a distinct planar dimensional 
orientation. They found a tendency towards a weak linear orientation of the 
c-axes of the pyroxenes in this plane and from these observations concluded 
that the fabric of this meteorite suggests it was formed by the settling of 
plagioclase and pyroxene crystals to the floor of a magma chamber where they 
accumulated in nearly horizontal layers. Consequently in the Moore County 
eucrite the plane of layering of the crystals approximates a plane perpendicular 
to the radius of the primary body at the site of crystallization of the meteorite. 
If it is assumed that the magnetic field in the upper mantle of the primary 
body was dipolar, like the Earth’s field, then the angle between this plane and 
the direction of magnetization in the meteorite will be the magnetic inclination 
of the field at this place and can be used to calculate the magnetic latitude of 
the site of crystallization of the meteorite. 

The direction of magnetization of a slice cut from the Moore County 
eucrite was determined by a well established procedure (Collinson, et al., 1957) 
and then the same slice was made into a thin section so that the orientation of 
the place of crystal layering could be determined relative to the direction of 
magnetization. The orientation of the c-axes of 34 pyroxenes was détermined 


on the universal stage and the results are plotted stereographically in figure 1. 
The c-axes of the pyroxenes lie in a plane at 90° to the plane of the thin sec- 
percent are {100} sections and 11 percent {001} sections. However the most 
important observation is thate the poles of almost all the sections lie in the 
plane of layering LL’ given by the pyroxenes, confirming a planar dimensional 


John Lovering 


axes of 34 pyroxenes relative to the plane of 
of magnetization M to the plane of layering LI 


orientation of the pyroxene and plagioclase crystals in the Moore County 


eucrite, A minor point of interest is the observation that the preponderance of 


{O10} sections in the thin section studied indicates a tendency towards a linea- 
tion of the plagioclase crystals. Hess and Henderson (1949) observed a similat 
weak lineation of the longest axes (i.e. the c-axes) of the pyroxenes, They also 
tion, This plane of layering would appear to correspond to that found by Hess 


ind Hendersor 1949) although our section seems to have been cut at about 


90° to the plane of their section 


In figure 2. the orientation of 35 plagioclase crystals in the thin section 


has been plotted, from which it ippears that 66 percent are fO1O! sections, 23 


found that the b- and ixial directions in the plagioclase crystals were elon- 


ited so that t of the crystals were lying with their {100} faces in the plane 


ol layering tions found in the present study are in accord with this 


conciusiol 


O10 
lOO 
OO! 


Poles of the {010}, {100}. and {001} plagioclase faces plotted on the plane 
ction. (LL’ is the plar 


plane of layering of the e-axes of the pyroxene defined in 


74 
ait 
ft 
: 
Fig. 1. Stereo projection of 
wil P is ) 
be 
3 9 
oe e 
° 
Fig. 2 
of the thin s 
figure 1) 


The Magnetic Field in a Primary Meteorite Body 275 


The direction of magnetization (M) in the thin section has been plotted 
in figure 1 and is inclined at an angle of about 20° to the plane of crystal 
layering in the section, giving a “geomagnetic” latitude of either 10°N or 10°S 
for the site of crystallization of the Moore County eucrite in the primary body, 


CONCLUSION 

The remanent magnetization observed in the Moore County eucrite sug- 
gests that while the meteorite formed part of the outer silicate mantle of the 
primary meteorite body it cooled through the Curie temperature of its most 
important magnetic constituent magnetite (560°C) under the influence of a 
magnetic field. Furthermore the Moore County eucrite would have cooled be- 
low this temperature (i.e. 560°C) while the metal core of the primary body 
was still at least partly fluid and capable of generating a magnetic field, From 
these conclusions it can be estimated that the temperature within the core at 
this time must have been at least 1700°C. 

The site of crystallization of the Moore County eucrite within the outer 
silicate mantle of the primary body can be estimated in two ways. Firstly, the 
existence of tridymite in the meteorite suggests that it probably crystallized 
within about 10 km of the surface of the body, Secondly, the inclination of the 
direction of magnetization in the eucrite to its plane of crystal layering (i.e. the 
horizontal plane at the site) suggests that the meteorite formed close to either 
the LO°N or the 10°S magnetic latitude of the primary body. 
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PETROFABRIC STUDY OF CONESTOGA LIMESTONE 
FROM HANOVER, PENNSYLVANIA* 
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ABSTRACT e stu applies established techniques of calcite petrofabric analysis toa 
ow ide met phic 11 textured, carbonate rock. Prominent features of the fabric 


osterystalline late phase orthorhombic symmetry, and 
fold axes. The fabric results from differently 


t the axes of which rotated through the rock 
causing late phase deformation are in the ac plane 


normal to slaty cleavage 


veme 


INTRODUCTION 


Experimental deformation of single crystals and aggregates of crystals of 
calcite has provided much data about the mechanical behavior of calcite under 
stress (Turner and Clvih. 1951: Borg and Turner. 1953). New methods of 
calcite petrofabric analysis and interpretation of fabric of naturally deformed 
marbles have resulted (Turner. 1953: MelIntyre and Turner, 1953). The 
present work was done to: 1. extend these methods of analysis and interpreta- 
tion to rocks of a different area that are of lower metamorphic grade than those 
previously studied; 2. describe and interpret the fabric of the Conestoga lime- 
stone. The Conestoga limestone was particularly suitable for this study because: 
1. megascopic structures such as slaty cleavage, and bedding are well-developed 
at the outcrop and provide optimum conditions for correlating field and micro- 
scope data (fig. 1); 2. the rock is a coarse-grained calcarenite with calcite 
grains large enough for easy optical measurement (fig. 2). Most other lime- 
stones sampled in the Great Valley and western Piedmont proved too fine- 
grained for universal stage work. 
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LOCATION AND DESCRIPTION OF OUTCROP 


The outcrop of Conestoga limestone sampled is 1.1 miles N 52 E of Han- 
over, Pennsylvania in a small abandoned quarry. It is the only outcrop in an 
8 square mile area mapped as Conestoga limestone by Weaver (1955), The 
rock is a limestone conglomerate (boulders to 1 foot in diameter) with sandy 
beds and matrix of coarse-grained quartz and calcite. Other writers (Wise, 
1953) described conglomerate near the base of the Conestoga limestone, This 
is the probable stratigraphic position of this conglomerate. 
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Fig. 1. Outcrop showing slaty cleavage (parallel to pencil) and bedding . 


Southeast dipping bedding (S,) is intersected by more steeply dipping 
slaty cleavage (S.) (fig. 1). Pebbles and boulders of limestone are flattened in 
the slaty cleavage and elongated down the dip of the cleavage in the a fabric 
direction. Phyllites in the area contain a mineral assemblage characteristic of 
the green schist facies of metamorphism. 


PETROGRAPHIC DESCRIPTION 

Specimens for this study were selected from the calcarenite matrix and 
beds. Thin sections were cut in the fabric ac plane, perpendicular to intersec- 
tions of bedding and slaty cleavage. In this plane all grains are elongated paral- 
lel to the trace of slaty cleavage, the ab fabric plane, and flattened normal to 
the cleavage (fig. 2). Long dimensions of grains measured 5 to 2.5 mm. 
Quartz was studied in both ac and be cuts. The mean ratio of long axes 
(fabric a) to short axes (fabric c) in 50 grains was 3:1 for both calcite and 
quartz. In a be cut the ratio of long axes (fabric 6) to short axes (fabric c) 
was 1.9:1 and 2.2:1 for calcite and quartz respectively. 

Essential components of the rock are quartz grains (38 percent), grains of 
fine to coarse crystalline calcite (40 percent), and aggregate grains of micro- 
crystalline calcite (22 percent). Quartz is slightly strained, replaced at all 


boundaries by crystalline calcite, and contains deformation lamellae in 50 per- 


cent of the grains observed. Crystalline calcite is clear, fractured into slivers 
parallel to slaty cleavage and cut by one or more sets of deformation lamellae 
(see below). This material was originally clastic, sand-sized crystalline calcite 
(e.g. fossil fragments) and clear calcite cement, Aggregate grains of micro- 
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Photomicrograph wing quartz, sand-sized crystalline calcite grains with 
deformation twins, and sand-siz regate grains of micro-crystalline caleite (NW and 


SW corners). Slaty clea s oriented NI SW. White line represents 0.5 mm. 


crystalline (.001 to .1 ite are sand-sized clastic fragments of very 
fine-grained lime mud 


THOD OF PETROFABRIC STUDY 


Measurement of orientation data in thin sections from two specimens (S-3 


ind S-18) was done on a Leitz 4-axis Universal Stage. Data were plotted on a 
| 
Schmidt-l imbert equal area projection net ol 20 cm diameter. 


Using procedures described by Turner (1949, p. 595-602) the following 
data were measured in each lcite grain: 


Orientation of the optic axis crystallographic axis |OOOL], 


?. Orientation of poles of all sets of visible deformation lamellae. In these 


specimens they were all found to be lamellae parallel to {0112} e. 
following plots were made from the measured data: 
1. Orientation of direction and sense of gliding in each {0112} lamella. 
Glide line or (O112) : 
Orientation ge [eze’| between any two {0112} lamellae in 


one grain 
Axes of compression ind tension most favorably oriented to have 
formed each {O112! twin. These directions are inclined at 45° to the 


slide line :r |. oriented so that gliding along the twin plane e would 
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occur in the sense required to produce twinning. (Turner, 1953, p. 
282. ) 

It is important to differentiate {0112} lamellae into descriptive classes 
that are comparable with lamellae produced in experimental work (Turner, 
1953; McIntyre and Turner, 1953; Borg and Turner, 1953, p. 1344-1346). 
When viewed in vertical position some lamellae of {0112} orientation appear 
as fine. dark lines that are not broad enough at 300 X magnification to show 
symmetrical differences in extinction with the parent grain. These are termed 
nontwinned lamellae. Other lamellae are clearly recognized as twins and are 
called twinned lamellae. Calcite grains studied have: no lamellae of any kind, 
one set of either nontwinned or twinned lamellae, 2 sets of either twinned 
lamellae or nontwinned lamellae. one set of twinned and one set of non-twinned 
lamellae. Table 1 shows the frequency of occurrence of the above classes, Al- 
though thin sections of only one orientation were studied, virtually all lamellae 
were seen because c-axes lie near the plane of the section. 


TABLE | 
Percent Occurrence of Twinned or Nontwinned {0112} Lamellae 
in Two Specimens Conestoga Limestone 
rype and number of sets of lamellae 
in each grain Spec. S-3 

No lamellae of any kind 
1 set, nontwinned lamellae 
2 sets. nontwinned lamellae 
1 set, twinned lamellae 30.! 


2? sets, twinned lamellae 


5. 
1 set twinned, 1 set nontwinned lamellae 2 12. 


100° 100.0% 


DESCRIPTION OF DIAGRAMS 

Kach writer analyzed one specimen independently and figures 3-12 show 
stereographic plots of data. Orientation of the plane of thin section is shown by 
the central strike and dip symbol in each diagram, Field orientation of the 
northwest end of the strike line is given by azimuth where 0 N, 90 E, 
180 S. 270 W. Angles of dip from 1-90° represent planes dipping to- 
ward the observer and angles over 90° indicate planes dipping away from the 
observer (Cloos and Hietanen, 1941, p. 37). Traces of bedding (S,) and slaty 
cleavage (S.) and positions of a. b, and c fabric axes are shown. 

Both sper imens show a well developed ac girdle of calcite c-crystallographic 
axes (figs. 3 and 4) and of £0112) deformation lamellae (figs. 5 and 6). S- 


C-axis maxima occur near the trace of the bedding and also at a large angle to 
the traces of bedding and cleavage. In both specimens, deformation lamellae 
maxima occur in pairs approximately bisected by the normal to slaty cleavage 
(S.). In S-18 an additional lamellae maximum is located in the bedding (S,). 
Figures 7 and 8 show central maxima for [e:e’] edges. That these maxima are 
close to the fabric 6 axis is at least partially a consequence of the ac girdle of 


Sets 
5-18 
shows better concentrations of both c-axes and deformation lamellae than S-3, 
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® 
Specenen S-5 Figure 3 Specimen S-16 Figure 4 

OO ¢c-anes, % Tl ¢-ones, contours- 96,7,42,14 % 


Specimen S-5 Specimen S-i6 Figure 6 
33 onto ? 84 lomelice, contours 72,47,24,12% 
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Specimen $-3 Figure 7 Specimen Figure 6 
37 e:e’ edges, contours- 108,6,54,26 % 18 @ e’ edges, contours- 22,16 11,5 % 


Figs. 3-8. Calcite c-axes, deformation lamellae and e:e’ edges, Note: Figures 3-14 


ire lower hemisphere projections of fabric data, Contours represent percent concentration 


of data in 1 ireas 
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c-axes and deformation lamellae. The more concentrated [e:e’] edge maximum 
of S-18 (fig. 8) is correlated with a narrow ac band of lamellae poles (fig. 6) 
while the more disperse [e:e’] edge diagram of S-3 (fig. 7) is correlated with 
a wide ac band of lamellae poles (fig. 5). 

Figure 9 is a plot of 9 grains in which two twinned sets of deformation 
lamellae were observed in specimen S-3. There is a strong grouping of c-axes 
of these grains in the zone perpendicular to the direction of applied compres- 


Specimen S-3 Figure 10 
4 


iv 


© Right-letere! glide lines 
© Left-ieterat glide Wines,//A 


Compression exes 
© Tension axes 


Specimen S-3 Figure Specimen S- 16 Figure 12 


Figs. 9-12. Caleite glide lines, compression axes and tension axes. 


sion shown in figure 11. Following Turner (1953, fig. 7) a plot of the bisector 
of the acute angle between glide lines [e:r] of adjacent twinned lamellae is also 
shown on figure 9, These bisectors lie close to the mean direction of applied 
compression shown in figure 11, Glide lines and sense of movement along them 
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404 


are plotted in figure 9, where it is apparent that the mean compression direc- 
tion (fig. 11) and grouping of bisectors separate zones of right lateral and left 
lateral twin gliding. Glide line orientation and relative slip direction along the 
elide line for all {0112! lamellae (both twinned and nontwinned) are shown 
in figure 10. Relative slip along glide lines has been classified as right lateral 
or left lateral with respect to the plane of the section, and the percent concen- 
tration of right lateral and left lateral gliding in two opposite quadrants 
hounded by diameters passing through mean compression and tension positions 
(fig. 11) has been computed. In quadrants Il and IV, 79 percent of the right 
lateral and 31 percent of the left lateral glide lines occur. while in quadrants | 
and Il, 69 percent of the left lateral and 21 percent of the right lateral glide 
lines occur (fig. 10). 

The stress that would be most effective in causing twin gliding on each of 
the twinned lamellae in the two spec imens is plotted in figures 1] and 12 as 
mutually perpendicular compression and tension axes, according to the method 
of Turner (1953, p. 282). Both diagrams show a broad zone of compression 
axes in a partial ae girdle symmetrically distributed about the normals to 
bedding and cleavage. Tension axes are most concentrated near the periphery 
in a narrower zone cul hy the traces of cleavage and bedding. In figure 11 
compression axes approach the center of the diagram from both peripheral 
concentrations suggesting a weak be girdle. Maximum areas of concentration 


ol compression and tension axes are marked C and T on both diagrams, 


INTERPRETATION OF DIAGRAMS 

Both indirect componental movements (recrystallization under stress) and 
direct componental movements (intragranular gliding. rupture and shear of 
grains, and rotation of grains) were probably effective in shaping the present 
fabric. The efficacy of each throughout the early stages of deformation can only 
he estimated: the role of posterystalline. direct componental movements neat 
the end of deformation is shown by petrofabric analysis, Calcite lattice orienta- 
tion developed before the postcrystalline, direct componental movements which 
are recorded in the diagrams of deformation lamellae and glide lines, Lattice 
orientation may be the result of early sedimentary or diagenetic orientation. 
recrystallization under stress, or early direct componental movements, evidences 
of which have been destroyed by post-lec tonic rec rystallization. The possibility 
that the major part of the calcite lattice orientation is the result of sedimentary 
or diagenetic processes is ruled out by the marked relation between fabric axes 
and c-cryst ulographic axis orientation. Sedimentary or early diagenetic orien- 
tation may be responsible for the c-axis maximum in both figures 3 and 4 near 
the trace of bedding (S,). Sander (1951. p. 6. 154) states that paradiagenetic 
calcite grows with c-axis parallel to bedding in Alpine Trassic limestones, The 
girdle of c-axes in the bedding produced by such growth could have been 
destroyed by later deformation except where the primary bedding girdle inter- 
sects the ac fabric girdle. thus giving rise to the maximum in the trace of 
heddine 

In strong experimental deformation the c-axes of calcite are brought into 
an apparently stable orientation at an angele of 10°-30° to the compression 
axis, (Turner, et al., 1956, p. 1284). This is the result of direct componental 
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movements without recrystallization, It has also been shown that syntectonic 
experimental recrystallization at 400° to 600° C. tends to align the c-axes of 
calcite parallel to the compression axis. (Griggs, et al., 1958, p. 48). 

In S-3 and S-18 compression axes sweeping through the rock perpendicu- 
lar to fabric 6 and passing through a number of different orientations at differ- 
ent times could explain the distribution of c-axes of calcite in an ac girdle 
(figs. 3 and 4). 

Alternatively, radial compression perpendicular to fabric b may have 
caused the ac girdle of c-axes. MacDonald (1957) has predicted on thermo- 
dynamic grounds that radial compression should form a girdle of calcite axes 
in the plane of compression. However, the presence of maxima within the ac 
girdle seems to indicate a number of different episodes of compression in 
different orientations as suggested above. The effect of recrystallization under 
stress in producing preferred orientation of calcite lattices is not known. 

Posterystalline, late stage. deformation of the rock has produced deforma- 
tion lamellae that may be interpreted dynamically, That this late deformation 
has had little effect upon pre-existing orientation of c-axes of calcite is shown 
in figure 9, C-axes of 9 grains showing 2 equally developed sets of twinned 
lamellae lie in a belt nearly perpendicular to the mean compression axis shown 
in figure 11. Obviously the original lattice orientation was the main factor 
controlling intensity of twin gliding. Other data shows that c-axes of grains 
with no deformation lamellae of any kind are all located in the direction of 
compression, a position where twin gliding on any of the potential {0112} twin 
glide planes would be impossible. 

Although Turner and Ch’ih (1951, p. 902) have shown that concentra- 
tions of [e:e’] edges parallel axes of extension in experimentally deformed 
marbles, central maxima of [e:e’] edges in figures 7 and 8 do not necessarily 


warrant this interpretation. A peripheral girdle of c-axes would be expected to 


sive rise to a central maximum of [e:e’] edges since ¢ A e:e’ 70°. Close cor- 
relation between width of ac girdles of deformation lamellae poles and c-axes 
and spread of b maxima of [e:e’] edges when comparing figures 3, 5 and 7 
and 4. 6 and 8 suggests that the degree of preferre d calcite lattice orientation 
of the earlier fabric is largely responsible for [e:e’] edge concentrations, 

Concentration of compression axes in a broad zone approximately per- 
pendicular to cleavage (S.) (figs. 11 and 12) agrees very well with observed 
flattening of grains, pebbles. and boulders in outcrop, Grouping of right lateral 
and left lateral glide lines show (fig. 10) the sense of movement that would be 
expected for the compression and tension axes plotted in figure 11, It is diffi- 
cult to illustrate in a 2-dimensional figure the system of 3-dimensional move- 
ment along glide lines, but since glide lines most commonly plot near the 
periphery of the diagram the writers believe that most twin gliding had a 
large component of movement in the ac fabric plane. The distribution of com- 
pression and tension axes (fig. 11). the distribution of right lateral and left 
lateral glide lines (fig. 10) and the distribution of double sets of stronely 
twinned lamellae (fig. 9) have an orthorhombic symmetry which is superim- 
posed upon the dominant monoclinic symmetry of c-axes and all deformation 
lamellae. 
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SYMMETRY AND INTERPRETATION OF FABRI( 


Sander’s (1950. part I. p 


produced the fabric has been confirmed by experimental deformation of 
marble, metal aggregates and clay aggregates (Turner, 1957, p. 11-16). A 
monoclinic symmetry of c-axes and all deformation lamellae is shown in figures 
». 4, 5 and 6. The symmetry plane is ac and 6 is the axis of symmetry, The 
symmetry. strong b-axes, and accompanying ac girdles mark both specimens as 
h B tectonites (Turner and Verhoogen, 1951, p. 548) except that no mega- 
scopic lineation 6 is recognized. B-tectonites are classically interpreted to result 
from either deformation on s-surfaces of shear intersecting in B or flexure fold- 
ing and rotation about B. in which case there are no intersecting (hO1) s- 
planes 

In the specimens studied, monoclinic symmetry of c-axes is the product of 
early stages of rock formation and deformation. As orientation of deformation 
lamellae hinges upon c-axis orientation the monoclinic symmetry of the plot of 
all deformation lamellae is indirectly also a product of early stages. although 
the actual formation of lamellae occurred during late phase, postcrystalline 
direct componental deformation. Orthorhombic symmetry is apparent in dia- 
grams of selected twinned deformation lamellae (fig. 9), directions of slip 
along all glide lines (fig. 10) and compression and tension axes of twinned 
deformation lamellae (figs. 11 and 12), all results of late phase, posterystalline. 
direct componental deformation, It is concluded that the monoclinic symmetry 
of the fabric is the product of successive phases of orthorhombic symmetry and 
movement that swept through the rock to different orientations during its his- 
tory of deformation. The last orthorhombic phase is preserved in unrecrystal- 
lized cle formation lamellae whe reas evidence of earlier phases has been ob- 


literated by post-tectonic recrystallization, 


The above discussion has dealt with distortion in the ae fabric plane only 
and be thin sections would have to be studied to discover distortion parallel to 
fabric b. For the following reasons. however, it is believed that distortion paral- 
lel to b has not been great: 

1. Macroscopic elongation of grains parallel to a, down the dip of the 

cleavage is visible in outcrop and hand specimen, Shortening parallel 


to b is considerably smaller. 


Most glide lines (fig. 10) fall near the ac plane and movement along 


them has a large component in the ac plane. 

Compression axes in at least specimen S-3 (fig. 11) approach the 
center of the diagram in a discontinuous be girdle, ruling out the pos- 
sibility of extension in 6, However. as suggested above, b could be a 


secondary axis of compression giving rise to the observed shortening 


in b 


Cloos (1947, p. 888) has shown in the South Mountain fold to the west 
that deformation in 6 is variable from place to place but generally small when 
compared to deformation in ac 


28 
ee. 31-83) basic assumption that symmetry of 
fabric of a deformed rock reflects the symmetry of movements which have 
: 
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DESCRIPTION OF QUARTZ DIAGRAMS (figs. 13 and 14) 

Quartz in specimen S-3 was studied in two mutually perpendicular thin 
sections parallel to ac and be. C-axes and available deformation lamellae were 
measured in one hundred grains in each slide and data from the be thin sec- 
tion were rotated into the ac plane for presentation. While in calcite virtually 


Figs. 13-14. Quartz c-axes and deformation lamellae. 


all deformation lamellae were visible in the ac thin section, such is not the case 
for the less strongly oriented quartz. By studying two perpendicular thin sec- 
tions all possible quartz deformation lamellae could be observed. 

Quartz c-axes show no definable pattern of preferred orientation and it is 
suggested that the cushioning action of surrounding calcite has prevented de- 
velopment of clear orientation. Cloos and Hietanen (1941, p. 65 and diagrams 
D-37 and D-38) demonstrate that quartz surrounded by other minerals com- 
monly shows a poor orientation, 

Deformation lamellae in quartz are located in two well defined areas of 
concentration corresponding to maxima for calcite deformation lamellae, Areas 
of minimum concentration for calcite and quartz deformation lamellae also 
correspond. In the experience of the writers it is common to find better concen- 


trations of quartz deformation lamellae than of c-axes. No fixed crystallographic 


relation of c-axes and deformation lamellae was found; the angle between c- 
axes and the perpendicular to deformation lamellae varies from 0-90° with a 
maximum between 10-40°. Interpretation of these relationships will have to 
await experimental work on quartz deformation, Thus far quartz deformation 
lamellae have not been produced in laboratory experiments. 


SUMMARY AND CONCLUSIONS 
Petrofabric analysis of calcite in a low-grade metamorphic rock with well- 
preserved sedimentary structures and textures has revealed surprisingly good 
preferred orientation of c-axes, deformation lamellae, slip along glide lines, 
and compression and tension axes. Petrofabric techniques of Turner (1953) 


& Gy \ \ \ \ 
Jo-' 
10 rotetea trom be cut ? Deformation Lamelice quartz. 73 trom ac cut, 34 
rotated from ut 
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ind others are definitely applicable to rocks of this metamorphic grade and of- 
ler excellent opportunities for the study of coarse-grained carbonate rocks of 
the Great Valley and western Piedmont. The monoclinic symmetry of c-axes 


ind deformation lamellae in the Conestoga limestone is a cumulative effect of 


many orthorhombic episodes of movement, differently oriented with respect to 


present fabric a and ¢ but always perpendicular to b. The last phase of ortho- 


rhombic. postcrystalline, direct componental deformation is recorded in de- 


formation lamellae analyzed in the present petrofabric study. 
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THE SYSTEM 
AT TEMPERATURES UP TO 200° 
A. J. ELLIS 
Dominion Laboratory, D.S.I.R., Wellington, New Zealand 


ABSTRACT. The equilibrium conditions in the system above were obtained by a dif- 
ferential vapor pressure method and the results up to 200° are reported. 

From measurements of the differences in vapor pressure between pure water and a 
salt solution it is possible in principle to obtain activity coefficients for the ions in solu- 
tion, The vapor pressures of sodium chloride and sodium carbonate solutions up to 200 
are given, but it was not possible to obtain results of sufficient accuracy to derive activity 
coethcients. 


INTRODUCTION 

Information on the system NasCO,—NaHCO,—CO,.—H.0O is required in 
understanding many problems which involve carbonate deposition in nature. 

Results giving the equilibrium conditions were reported by Walker, Bray, 
and Johnston (1927) at 25° and 37°, and more recently by Mai and Babb 
(1955) at temperatures up to 65°C. Both workers used a wide range of con- 
centrations, and the agreement between the two sets of results is satisfactory. 

A sensitive method of measuring small differences in the vapor pressures 
over two liquids is often required in the investigation of two phase aqueous 
solutions at high temperatures (e.g. equilibria involving a second gas, accurate 
determination of the vapor pressure of salt solutions). Results from the present 
investigation showed that the mercury-in-glass manometer could be used to 
record small pressure differences up to at least 250°. The limiting factor in 
measuring small differences in water vapor pressure is the minimum tempera- 
ture difference that can be maintained between the two containers. The be- 
havior of a suitable differential pressure apparatus was tested on the carbonate 
system and in determining the vapor pressure of salt solutions. 


THEORY 
In this paper carbonic acid refers to the molecule H.CO, and not to the 
mixture (H.O + CO.) often called by this name, At 25°, 0.259 percent of the 
total carbon dioxide in a pure water solution consists of H.CO, (Wissbrun, 
French, and Patterson, 1954). 
For the liquid-vapor equilibrium of carbon dioxide the following equation 
holds. 
(CO. ) (1) 


where B is an inverse Henry’s Law coefficient, (CO.) is the activity of molecu- 


lar carbon dioxide in solution, and f,,, its fugacity in the vapor. In this paper 


rounded brackets indicate activities of species, and square brackets concentra- 
tions. 
For the first and second ionization constants of carbonic acid at constant 
temperature and pressure the following relationships exist. 
(H+) (HCO,”) /(H.CO; ) kK, (2) 
(H+) (CO,”) /(HCO,’) kK. (3) 


The hydration of carbon dioxide in solution to carbonic acid is given by 


9707 
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the equilibrium constant 
K. 
where fy. is the fugacity of water vapor. For the ionic reaction 
2HCO,’( soln.) ~CO,” soln.) + gas) + gas) 
the equilibrium constant K is equivalent to 
K (CO;”). (HCO,’)? 
By substitution from equation (1) to (4) it can be shown that 
K K./K,. B. K, 

The quantities determined directly from the experiments are the concen- 
trations of sodium carbonate and bicarbonate in solution and the partial pres- 
sures of water and carbon dioxide. 

Let poo. p NaHCO, |? k, 
(HCO,’) = y, [NaHCO,}, 
and y2 [Na-CO,] . 
y: and ys are the respective ion activity coefficients for the HCO,’ and CO,” 


ions, k may then be related to K. 
k K. 


kK. 
a, and a, are the fugacity coefficients of CO. and H.O in the vapor phase. 
The apparent acid dissociation constant of the mixture (CO, + H,O) is 
labelled K,’. 
Ky (H+) (HCO,’)/(H.CO, + CO.) 
(H*+) (HCO,’)/(CO.) (1 + 
(1 + K,. fino) 
The error involved in the last approximation is equal to the percentage of 
total carbon dioxide present in solution as carbonic acid. This is unlikely to be 
greater than | percent in the temperature range examined. 


THE PRESENT EXPERIMENTAL METHOD 
The present study of the system Na.CO,—NaHCO,—CO.—H.O, was 
made by a differential pressure technique. Figure 1 is a diagrammatic outline 
of the apparatus. In a 3 inch stainless steel block two holes were drilled; the 
large one had a volume of approximately 105 ce and the other one held ap- 
proximately 5 cc. In the experiments the small chamber contained pure water, 


and the large one a carbonate solution of known initial composition, The dif- 


ference in pressure, as indicated by a thick-walled glass manometer, was equal 
to the partial pressure of carbon dioxide above the carbonate solution less the 
lowering of water vapor pressure by the salts in solution, The glass manometer 
was shielded by heavy metal gauze in case of an explosion. 

The technique was as follows: the high pressure lines (capillary stainless 
steel tubing) were removed from the top of the glass-metal connection and 
connected to dead-end pieces, The bombs and all the lines at this stage were 
completely dry. Known quantities of sodium bicarbonate and carbonate were 
placed in the large chamber, and both heads placed on the bomb. 

With all valves open the system was quickly evacuated to better than 1] 


at Temperatures up to 200° 


DIFFERENTIAL 
CAPILLARY VAPOUR 
TUBING PRESSURE 


APPARATUS 
FIG. | 


Vi 


GLASS-METAL 
CONNECTION 


GLASS 
TEFLON | MANOMETER 


Fig. 1. Apparatus for measuring small differences between the vapor pressure of two 
solutions. 


mm vacuum by a rotary oil pump. This was found not to affect appreciably 
the dry NaHCO,. The evacuation was through valve 4 which was also con- 
nected to a source of degassed water through a T piece. 

The degassed water was prepared by boiling water under reflux for 2-3 
hours under a partial vacuum in an apparatus of the type described by Law 
(1949). With valve 2 closed, approximately 70 cc of water was let into the 
large chamber and valve 1 closed. Valve 2 was temporarily opened to admit 
a few cc of water into the smail chamber. At this stage all the lines were filled 
with degassed water at atmospheric pressure. The manometer was always filled 
with water from the mercury surface up to the top of the glass-metal connec- 
tion. With the latter overflowing with water the lines were again connected to 
it with valves 3 and 4 closed. 

The bomb was heated in an 18 inch diameter 10 gallon bath of oil which 
was surrounded by a foot thickness packing of diatomaceous earth, A thermo- 
stat controlled the temperature to about + 0.1° and the bath was vigorously 
stirred. The fact that the two chambers were in the one massive piece of metal 
damped out temperature variations in the bath. With pure water in both 
cavities the differences in pressure registered at 200° were of the order of 
+ 0.3 em. As at 200° the change of vapor pressure is about 24 cm for a 
degree temperature difference, the difference in temperature between the two 
chambers was less than 0.01°. 
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After the apparatus h id been at the temperature of observation for about 
in hour, valves 1 and 2 were opened and the values of Ap read on the manom- 
eter until they remained constant. This usually required less than an hour. 
With a little experience the manometer could be pre-set to about the value of 
Ap expected so that there was no appreciable transfer of water into either the 
solution or the small chamber when the mercury moved. 

The weight of solution in the bomb was measured after the experiments 
ind » concentrations estimated. In the calculations, corrections were made 

water transferred by changes in the mercury level. the water in the 
vapor space, the depression of the water vapor pressure by the salt in solution, 
the carbonate in the bicarbonate used. and the change in the original car- 
bonate-bicarh itio by the reaction to form carbon dioxide in the vapor 
and liquid 

The only correction whi h was very appre iable was the depression of 
the water vapor pressure by the salts in solution, Values for this correction 
were taken from experiments reported below with the assumption that the 
osmotic effect of NaCl equalled that of NaHCO,. 

The lines from the bomb were arranged symmetrically about the manom- 
eter and the heads of the two chambers were at the same level. The values of 
Ap were corrected to pressures in centimetres of mercury by taking the effective 


density of the manometer liquid as 12.55 @ ‘ce. 


RESULTS 


s were completed at ionic streneths equal to 


Three series of ¢ \periment 


0.5. and 1.0. A typical experimental result is calculated in detail. 


| 162°C 

| NaHCO, used contained 0.0034 moles 

NaCO, mole NaHCO 
0.01063 moles 
0.0385 moles 
67.9 
187.7 cm Hg 
11.6 em Hg 
176.1 em He 
12.9 em Hg 
24.5 em Hg 
76.6 
28.8 


0.00048 


0.03755 
0.01111 
[NaHCO 0.553 moles/1000 water 
) 1.162 moles water 
[NaHCO 1.08 at. 
ke.mole 


Phe total quantity of carbon dioxide in the system was calculated from 
the partial pressure of the gas. the volumes of gas and liquid. and the solu- 
lity of carbon dioxide reported in the first paper of this series (Ellis, 1959). 
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The results from the experiments are given in table 1 and the values of 
log k are plotted against 1/T° K in figure 2. 

It was shown that K k. Y2-@)-Q2/ ¥ 
where K is the thermodynamic equilibrium constant for the reaction 

2HCO,’ (soln.) CO,” (soln.) + H.O(gas) + CO. (gas) 
At constant pressure and composition 
din K dln k Aln ( 
a7 oT ot 
AH' AH + Ah 

where AH® is the standard increment in heat content for the reaction and Ah 
incorporates the overall change in relative partial molal heat contents of ions 


T T 
SYSTEM  NegCO, - NeHCO, - 
COz = 


FIG. 2 


ABOVE 100°C | 
= 

x O-5 
° = 0.2 


1/T°K x 
3.2 30 28 £26 24 22 


rhe results obtained for k at various temperatures 
(k [NasCOs]. Py. /lNaHCOs]). 


in solution and gases in the vapor. Table 2 gives values of AH obtained from 
the experiments. 
2 
ree 25 50 100 150 200 
AH(u 1) 16400 16400 15800 14600 13000 cals./mole 


AH(x 0.2) 15900 15900 15500 14200 12500 cals./mole 
AH uo (Vapzn.) 10480 10230 9697 9074 8348 cals./mole 


AH* (xu 0.2) 5700 5670 5800 5130 $150 cals./mole 


AH (vapzn.) is the latent heat of vaporization for pure water (Keenan 
and Keyes, 1936), and assuming that this heat of vaporization is the same for 
a solution of ionic strength 0.2, AH* AH AH (vapzn.) is the 
experimental increment in heat content if the reference state of water is taken 
as the liquid rather than the gas. 

kor the reaction 

2NaHCO.(soln.) Na.CO,(soln.) + H.O(1) + CO.(¢) 
ACp* is very small at temperatures below 100° (AH* is constant within ex- 


> 
-| 
LOG k | 
“2 
ZA 
| 
J = 
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perimental error). At higher temperatures AH for vaporization of CO, from 
solution decreases rapidly (see first paper of this series) and this will be a 
major factor in the decrease of AH*. 


If as an approximation in the equation 


2 
K,.K..B 
«, are unity, and that K,’ K,.Ko.fio, 

Ky’.B. ye 
Approximate values of k’ are calculated in table 3 from known experimental 
data (K,’ up to 50°, Harned and Davis, 1943; Kk. up to 50°, Harned and 
Scholes, 1941; and up to 90°, Cuta and Strafelda, 1954), Values of Ky’ and K, 
are extrapolated to higher temperatures by means of equations of the type used 
by Pitzer (1937). Values of B were obtained from the solubility results re- 
ported earlier. 


, it is assumed that the fugacity coefficients a, and 


then k 


Mean ion activity coefficients Y-+ determined by experiment are related 


to the individual ion activity coefficients by the equation 
y'4.y"— (Harned and Owen, 1950, p. 10) 
where n and m are the number of positive and negative ions respectively pro- 
duced hy a salt. 
3 
B 
x 10? Puyo k 
(mol, (at.) k’ (u=1.0) (u=1.0) | 
3.43 0.0313 9.610 26x10" 
1.89 0.1217 8.4x 
1.18 1.00 11x10 6.310? 
1.14 4.70 1.4 10° 0.74 r 
27 15.3 0.16 4.3 27 


(NaHCO, soln.) 


Therefore 
(Na.CO, soln.) 


(NaCl soln. ) 
(NasCO, soln.) 


Values of y are reported by Harned and Owen (1950) p. 557, for NaCl 
solutions up to 100°, and for Na.CO, solutions up to 95° by Taylor (1955). 
For solutions in the series y 1.0, [NasCO,] and [NaHCO,] equal approxi- 
mately 0.1 m and 0.75 m respectively. Values of the ratio, 

* (NaCl soln.) y.* (NasCO, soln.) 


at these concentrations are given in table 4. 


TABLE 4 
te 25 50 75 125 
y> (0.75m NaCl) 0.666 0.665 0.643 0.600 
(0.1m NazCOs) 0.465 0.524 0.505 0.30 
Cale. Ratio y:°/4 1.95 1.4 1.3 4.8 
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The calculated values of y,°  y. in table 4 are intended only as a qualita- 
tive illustration and cannot be expected to equal the experimental values in 
table 3. The mean activity coefficients for 0.lm Na.CO, in a 0.75m NaHCO 
solution, will be different, and probably lower than y. (Na2CO,;) in a pure 
0.1m solution. 

However, the figures show that the rapid increase in y,"/y2 at higher tem- 
peratures is caused mainly by the decrease in the mean activity coefficient of 
sodium carbonate. This would indicate that this substance is only slightly 
dissociated in high te mperature water, Its decrease in solubility to almost zero 
at the critical temperature of water is therefore not unexpected. 

VAPOR PRESSURES OF NaCl and Na.CO, soLUTIONS 

The water vapor pressures over solutions of sodium chloride and sodium 
carbonate were determined up to 200° using the differential pressure apparatus. 

The procedure used was to dry and evacuate the bomb and lines, and 
with valve 1 closed. to allow degassed water to enter the small chamber. Valve 
2 and valve 4 were then closed. The solution to be studied was also degassed in 
the reflux apparatus and approximately 70 ce of degassed solution was allowed 
to enter the large chamber through valve 1, valve 3 and valve 4. All valves 
were then closed and the manometer connections made as before. The manom- 
eter was pre-set to the approximate level expected in the experiment. 

The concentration and weight of the solution was obtained after the ex- 
periment, and where appropriate, corrections were made as before. The practi- 
cal osmotic coefficient @ is defined by the equation 


950, p. 13). The activily of 


water a, 3 f where fy... and [no are the fugacities of water 


over the solution, and over pure water. If it is assumed that the fugacity ratio 
can be replaced by p p . Le., the fugacity coefficients are the same 
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Fig. 3. The difference between the vapor pressure of sodium chloride solutions, and 
pure water at various temperatures, 
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for water vapor over the solution and over pure water, then 
P 

t.m P 

where t is the number of ions formed by the salt, m its molality, pyjo the 
water vapor pressure over the solution, and p’,,. that over pure water, 


The results of the determinations are given in figures 3 and 4, Lines for 
various concentrations, obtained by giving ¢ a value of unity, are also shown. 

Harned and Owen (1950, p. 363) tabulate osmotic coefficients for sodium 
chloride solutions up to 100°. Table 5 gives some of their values for @ and the 
values obtained with sodium chloride in this work above 100°. 

Within experimental error the present results follow on smoothly from the 


previous values. 
5 


Values of ri) for NaCl Solutions 


100 120 140 160 180 200 


60 80 

0.9210 0.9178 0.9135 - 0.82 0.80 

0.9234 0.9195 0.9145 0.91 0.91 0.88 0.82 0.77 

0.9442 0.9402 0.9345 0.91 0.88 0.86 0.785 0.705 

It was planned originally to obtain mean activity coefficients for ions of 
various salts from the vapor pressure measurements by the standard process 
(Taylor, 1955) of applying an integrated form of the Gibbs-Duhem equation. 


a,, atm 55.5] 
d log a, —log m/m, 


5 (m) 
y+ (m,) a,;, at m, t.m 


where m, is a reference molality at low salt concentrations below which values 
of the activity coefficients can be extrapolated to infinite dilutions by the 


Debye-Hiickel theory, and a, is the activity of water. 
To derive useful values for y, (m) it would be necessary to obtain values 


+28 0.9 Na,CO, 
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of Ap for solutions down to about 0.1m with an accuracy of the order of 1 per- 
cent. It is unlikely that the present results are better than + 3mm in Ap over 
most of the range. Even at 200° this would amount to an error of almost 
10 percent in activity coefhicients at concentrations of 0.1m. 

The results are therefore left in terms of pressure differences and osmotic 
coeflicients, and as such should be useful in calculating water vapor pressures 
over hydrothermal solutions. 
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PHASE EQUILIBRIA IN THE SYSTEM 
MANGANESE OXIDE—SiO., IN AIR* 


ARNULF MUAN 


Dept. of Metallurgy, Pennsylvania State University, 
University Park, Pennsylvania 


ABSTRACT, Phase relations in the system manganese oxide—SiQz in air have been in- 
vestigated in the temperature range 990-1571°C, using the quenching technique. The data 
obtained permit the construction of a phase diagram showing stability relations in the 
liquidus temperature region as well as at subsolidus temperatures, For sake of simplicity 
the phase relations are illustrated in the form of a projection into a chosen join through 
the ternary system. The resulting diagram has the appearance of a binary system, with 
“binary invariant” situations as follows: a peritectic exists at 1230°C, with cubic MnsQ,, 
tephroite (2MnO-SiO.) and liquid coexisting in equilibrium with a gas phase of Oz par- 
tial pressure equal to 0.21 atm. At 1206°C tephroite, MnO- SiO, solid solution (rhodonite) , 
liquid and gas are present together in what appears as an eutectic situation, Tridymite 
(Si0.), MnO-SiOe(ss), liquid and gas exist together in equilibrium in an apparently 
peritectic situation at 1272°C, and cristobalite together with 2 liquids and gas are the 
equilibrium phases at approximately 1700°C, In the subsolidus region the following ap- 
parently invariant situations have been determined: at 1204°C cubic MnjQ,, tephroite and 
MnO-SiO2(ss) coexist in equilibrium with a gas phase of Oz partial pressure equal to 
0.21 atm. Tetragonal Mns0,, a phase probably analogous to the mineral braunite 
(MnO3(ss)), and MnO-SiO.(ss) are in equilibrium with this gas at 1168°C, while at 
1048°C) MnzOs(ss), MnO-SiOe(ss) and tridymite are the stable phases in air. In addi- 
tion, “invariant” situations also result from the tridymite-cristobalite inversion at approxi- 
mately 1470°C and the inversion from tetragonal to cubic MngQ, at approximately 1160°C. 


INTRODUCTION 


Phase equilibrium studies are generally recognized as being necessary 
first steps in any investigation aimed at improving our knowledge of hetero- 
veneous reactions. Accurate measurements of chemical properties in such 
systems can be interpreted intelligently only if the system is well defined with 
respect to the phases present. 

Oxide systems have been the center of interest of research workers in 
many different fields. Geochemists and petrologists have gathered a tremendous 
amount of phase equilibrium data which have enabled them to infer conditions 
governing the formation and existence of many minerals and rocks, Their 
studies have dealt mainly with combinations of oxides of noble gas type ions 
such as for instance SiO., Al.O,, CaO, MgO, Na.O, K.O, while data on most 
oxides of transition elements, such as for instance manganese oxides, are rela- 
tively scarce. This situation undoubtedly exists partly because the latter sys- 
tems are more difficult to study experimentally, but also because transition 


elements in their natural occurrences frequently are present as sulfides rather 


than as oxides. 

Oxides also play a very important role as one or more of the phases 
present in many high temperature industrial processes. Ceramists deal es- 
sentially with oxide materials, and their efforts have contributed significantly 
to our knowledge of oxide chemistry. It is natural that ceramists are interested 
primarily in the chemically very stable oxides of noble gas type ions, such as 


Contribution No. 58-10 from College of Mineral Industries, The Pennsylvania State Uni- 
versity, University Park, Pennsylvania. 


297 


298 {rnulf Muan—Phase Equilibria in the 


SiO.. Al.O., CaO. MeO, Na.O. K.O,. in other words the same oxides as are 
studied by geochemists. In the field of metallurgy, however, one is mainly 
concerned with the extraction of metals from oxide raw materials, and hence 
oxides of low stability are the center of interest. such as for instance iron 
oxides and manganese oxides. The meta!lurgist cannot, however, overlook the 
importance of systems made up from among the more stable oxides, which 
form the constituents of the slag phase and also are used in the refractory 
linings of the furnaces 

Phe most important oxide by far among those of transition elements is 
iron oxide, Systems containing iron oxide as a component have been studied 
and restudied by metallurgists. geochemists and ceramists alike over the last 
50 years. By comparison, very few studies have been made of systems where 
oxides of other transition elements are present. Among those of next im- 
portance both technologically as well as geologically is manganese oxide, 

While data on the system Mn—O are gradually accumulating in another 
phase of our research efforts (Hahn and Muan, in press), the present work 
represents a first step in our investigation of the system Mn—Si—9, Because 
manganese occurs in different states of oxidation, the O, partial pressure of 
the gas phase must be controlled carefully in any investigation of a system 
where manganese oxides are present. We have chosen to start our study of 
manganese oxide—SiO. equilibria by investigating the system in air, that is 
it a constant O, partial pressure of 0.21 atm. In choosing this approach we 
are actually working along the irregularly curved 0.21 atm, O. isobaric sec- 
tion through the system Mn——Si—O. Equilibria along this section may be 
represented in a simplified manner by projecting the points representing com- 
positions of condensed phases into a chosen join. We will present the results in 
terms of the components Mn.O, and SiQ,. The resulting diagram has the ap- 
pearance of a binary system. Apparently invariant situations in this system 
are actually intersections between the 0.21 atm, isobaric surface and univariant 
lines in the system Mn—Si—QO. Along such liquidus univariant lines two 
erystalline phases. liquid and gas coexist in equilibrium, whereas a subsolidus 
univariant situation is characterized by the coexistence in equilibrium of three 
crystalline phases and a gas phase. 


Previous Work 


Although manganese silicates have attracted considerable attention among 
geochemists as well as metallurgists, no systematic study of phase relations in 
the system manganese oxide—SiQO, in air has been reported in the literature. 

The one end member in the system, SiO., of course has been the subject 
of numerous phase investigations over the past few decades, starting with the 
now classic work in 1913 by Fenner. In spite of all efforts since that time, 
doubt still exists as to certain of the important characteristics of this substance, 
particularly with respect to the stability relations of the tridymite phase, An 
impression of the present status of our knowledge of this system can be ob- 
tained by reading recent papers by Sosman (1955), Flérke (1956), Eitel 
(1957), Hill and Roy (1958). 


The system Mn—O, unfortunately, has not been studied experimentally in 
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any great detail. Scattered equilibrium data and calorimetric data are avail- 
able in the literature, from which Coughlin (1954) has tabulated the standard 
free energies of formation of the various manganese oxides as a function of 
temperature, Most of the values listed in these tables were obtained by addi- 
tion and subtraction of large numbers, each of which is subject to considerable 
margins of error. Hence the stability relations of the manganese oxides de- 
picted by Coughlin’s tables are uncertain and should be considered only as 
rough approximations. Hahn and Muan (in press) have studied recently 
limited parts of the Mn—O system, particularly the equilibria between Mn,O, 
and Mn,O, and between Mn,O, and MnO. Their data show that Mn,Q, is the 
stable crystalline phase of manganese oxide in air above 877°C and up to the 
melting temperature of 1567°C. At temperatures slightly below 877°C Mn,O; 
is the stable phase in air. It is well known that Mn,O, exists in two different 
stable modifications, one high temperature cubic form with spinel structure 
and one low temperature tetragonal form, The exact temperature at which 
these two modifications coexist in equilibrium is still not certain. A discussion 
of the problem appears in a recent publication by Van Hook and Keith (1958). 
They indicate that the most likely inversion temperature is approximately 
1160°C, The oxide Mn.O,, likewise, has two established polymorphic forms, 
the a-Mn.O, with “C”-type cubic structure (Zachariasen, 1928; Pauling and 
Shappell, 1930) and the y-Mn.O, with distorted spinel structure (Dubois, 
1935). The latter modification, however, is unstable relative to the former and 
will not appear in our phase equilibrium diagram. 

Some of the previous experimental phase equilibrium studies of mixtures 
of manganese oxide and silica have suffered from insufficient control of the O. 
partial pressure of the gas phase in equilibrium with the condensed phases. 


This is probably one of the reasons why such a confusing variety of phase 


diagrams for the system “manganese oxide—silica” have appeared in the 
literature (Doerinckel, 1911; Glaser, 1926; Herty, 1930; White, Howat and 
Hay, 1934). A diagram for the system manganese oxide—SiO, under con- 
trolled, strongly reducing, conditions was published in 1958 by Glasser, who 
used CO./H. mixtures in order to define more precisely the O, partial pressure 
of the gas phase. 

\ large number of descriptions of naturally occurring manganese min- 
erals with a bearing on the present study are found in mineralogy textbooks 
(see for instance Dana’s System of Mineralogy, 1944) and numerous publica- 
tions. Papers of particular interest are those describing the mineralogy and 
geology of the Langban deposits in Sweden (Magnusson, 1924; Flink, 1926; 
Palache, 1929). 

Experimental Method 

General Procedure.—The phase relations were determined by the quench- 
ing method, Preheated mixtures of manganese oxide and silica were held in 
platinum containers in air at constant temperature until equilibrium was 
reached among crystalline, liquid and gas phases. The samples were then 
quenched rapidly to room temperature, Phases present were identified by 
x-ray and microscopic methods, and compositions of liquids at liquidus tem- 
peratures were determined by wet chemical analysis. 
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Materials.—Oxides of highest purity available commercially served as 
starting materials. The source of silica was c.p. silicic acid which was dehy- 
drated by heating at 1300°C for 12 hours. The source of manganese oxide was 
“Baker Analyzed” grade manganese dioxide, MnO.. Samples weighing 10 g 
were made up by mixing these materials in required proportions and me- 
chanically grinding the mixture. Each mixture was then heated slowly to 
melting either in air atmosphere in a globar furnace, or in a gas-air combus- 
tion furnace in case of higher melting mixtures. The heating of the starting 
materials had to be carried out with utmost caution in slow steps from low 
temperature up to the liquidus region because the successive decomposition of 
higher to lower manganese oxides took place very violently. Particularly 
troublesome in this respect was the temperature region in which the first 
liquids form. These preheated and quenched mixtures, after grinding to minus 
60 mesh, were used for the equilibration runs for determination of liquidus 
temperatures. 

The same mixtures were usually used for determination of boundary 
curves also in the subsolidus temperature region. As a check on attainment of 
true equilibrium a few runs in this region were made with two different start- 
ing materials prepared by prolonged heating of the mixtures, one at tem- 
peratures slightly above and one at a temperature slightly below that of the 
boundary curve. 

Furnaces and Temperature Control.—Vertical tube furnaces with wind- 
ings of platinum or 80 platinum-20 wt. % rhodium were used for most of the 


quenching experiments. A modified Roberts and Morey (1930) strip furnace, 
using strip resistors of a 60 platinum-40 wt. % rhodium alloy, was used for a 
few quench runs at temperatures above 1600°C. 


Temperatures in the vertical tube furnaces were measured before and 
after each run with a platinum 90—platinum-10 wt. % rhodium thermocouple 
frequently calibrated at melting points defined as follows: Au, 1063°C; 
CaMeSi.Og. 1391.5°C; CaSiO,, 1544°C, Literature values for temperatures on 
the Geophysical Laboratory scale are used without correction up to 1550°C, 
the upper limit of the original scale. Temperatures above 1550°C are adjusted 
to the 1948 International Scale using the correction data tabulated by Cor- 
ruccini (1949) and by Sosman (1952). Furnace temperature was controlled 
by a second thermocouple inserted close to the heating element and connected 
through compensating lead wires to a Celectray controller. The temperature 


MnsQ, and SiO, chosen as components. This part of the figure shows phases present as a 
function of temperature and manganese oxide to silica ratios of the mixtures. Heavy lines 
are boundary curves. These are dashed in regions where experimental data are less re- 
liable or lacking, Solid dots, solid triangle, open circles and crosses represent experi- 
mentally determined points on liquidus, 2 ag solidus and subsolidus boundary curves, 
respectively, based on data contained in table The lower, triangular diagram shows 
true compositions of liquid and crystalline + oft expressed in terms of the chosen com- 
ponents MnO, MnzOs and SiQs. Solid dots represent experimentally determined composi- 
tions of liquids at liquidus temperatures, and the heavy curve is the 0.21 atm. liquidus 
isobar drawn through these points. Light straight lines indicate approximate locations of 
conjugation lines connecting points representing compositions of liquids with points 
representing composition of crystalline phases with which the liquids are in equilibrium. 
Light dash curves are used to indicate how the MnsO,—SiO. join is swung up to the 
Mn-O;—SiOz join for projection purposes, as explained in the text. 
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gradient at the hot spot did not exceed 2°C over 1 inch at 1300°C, Tempera- 
tures on the strip furnace were measured with an optical pyrometer calibrated 
at melting points defined as follows: CaSiO,, 1544°C; 10 CaO 90 wt. % SiO., 
1708°C; Pt, 1769°C, 

Examination of Quenched Samples.—Phase identification in the liquidus 
temperature region was most easily accomplished with the microscope, using 
transmitted as well as reflected light techniques. The examination of polished 
sections in reflected light was particularly useful for distinguishing primary 
crystals from quench crystals in the manganese oxide rich part of the system. 
The identity of the phases was further confirmed by x-ray techniques, using Fe 
radiation with a Norelco spectrometer unit. The x-ray technique was also used 
for determination of compositions of solid solution crystalline phases in the 
subsolidus temperature region. 

In order to determine the compositions of liquids at liquidus temperatures 
it was necessary to carry out chemical analysis of quenched samples. After 
crushing of these samples, suitable amounts (~0.2 g) were weighed and trans- 
ferred to a platinum container. Excess standard ferrous ammonium sulfate 
solution was added together with 5 ml concentrated H.SO,, and the solution 
was boiled for a couple of minutes. After cooling, 10 ml of HF was added, a 
cover was put on the crucible and the solution boiled until the sample was 
completely dissolved, which usually required a period of approximately 10 
minutes. The crucible was plunged into a 600 ml beaker containing 250 ml 
water and 10 ml 1:1 H.SO,. After saturation with boric acid, the solution was 
titrated with 0.05 — n KMnQ,. 


RESULTS 


Results of quenching experiments are listed in table 1 and illustrated 
graphically in figure 1. This figure consists of 2 diagrams. The upper diagram 
with the appearance of a binary system shows phases present as a function of 
temperature. Because the compositions of all condensed phases cannot be 
represented in terms of the two chosen components, a supplementary diagram 
is presented in the lower half of the figure. This ternary diagram shows true 
compositions of condensed phases in equilibrium along boundary curves in 
the upper diagram. 

The most important features of the diagrams are as follows: liquidus 
temperatures are high at both ends of the system, 1567 and 1723°C being the 
melting temperatures in air of manganese oxide and cristobalite, respectively. 
In the silica rich end of the system, liquidus temperatures remain high 
(~1700°C) until a composition of 50 Wt.% manganese oxide is reached, 
corresponding to the maximum extent of liquid immiscibility, As the man- 
ganese oxide content increases above this value, liquidus temperatures fall off 
sharply. The MnO-SiO,(ss)-tridymite boundary curve is crossed by the 0.21 
atm. QO, isobar at 1272°C in what appears as a peritectic situation in the upper 
diagram of figure 1. The lowest liquidus temperature in air is 1206°C, with 
MnO-SiO.(ss) and tephroite coexisting with a liquid containing 66 Wt.% 
manganese oxide and 34 Wt.% SiO,. As the manganese oxide content increases 


above this level, liquidus temperatures start rising, and another “peritectic” 
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situation prevails at 1230°C, with tephroite and Mn,O, as the crystalline 
phases present in equilibrium with liquid. In the composition range between 
67 and 100 Wt.% Mn,O,, cubic Mn;O, is the primary crystalline phase. 

Tephroite is stable in air in the temperature interval 1204-1230°C. At the 
upper temperature limit of stability it decomposes incongruently to Mn;O, 
and liquid; at the lower temperature limit of stability it decomposes to Mn;O, 
and MnO-SiO.(ss). The phases Mn,O, and MnO-SiO.(ss) coexist in equilib- 
rium down to 1168°C in air, At this temperature Mn.O,(ss) becomes stable 
and the condensed phases present are either Mn,0, (ss) or Mn.O; (ss) 

MnO-SiO.(ss) at temperatures below 1168°C, depending on the SiQ, con- 
tent of the mixture. The MnO-SiO.(ss) phase exists in stable equilibrium in air 
down to 1048°C, Below this temperature it decomposes to Mn,O;(ss) and 
tridymite. Tetragonal Mn,O, exists as an equilibrium phase in air down to a 
minimum temperature of 877°C. Below this temperature pure Mn,Q, is oxi- 
dized to Mn.O,; this transformation temperature increases to a maximum of 
1168°C as the SiO. content of the mixture increases. 

A summary of “invariant” situations in the system manganese oxide— 
SiO. in air is presented in table 2, and powder x-ray diffraction data for 
manganese oxide containing crystalline phases are listed in table 3. 


DISCUSSION 


The most interesting feature of the phase relations in this system is the 
stable existence of a SiO.-containing solid solution phase with a—Mn,Q, struc- 
ture. The large increase in the transformation temperature of Mn,O, to Mn;O, 
as SiQ, is added to manganese oxide shows conclusively that SiO, enters the 
structure of the latter. The phase formed is termed Mn.Q,(ss) in this paper, 
leaving the question open as to the exact nature of the solid solution, The 
x-ray patterns show no appreciable shift in the d-spacings of the main lines of 


Mn.0, as SiO, content increases, but some weak extra lines appear, indicating 
some degree of ordered substitution rather than a completely random distribu- 
tion of the cations. 


The Mn.O,(ss) phase appears to be identical to the naturally occurring 
mineral braunite, which has been described and studied by a large number 
of investigators in the past. Structure determinations by x-ray diffraction 
studies have been reported by Aminoff (1931) and by Bystrém and Mason 
(1943). The structures as well as d-spacings are very similar to those for the 
pure Mn.O, end member, which has a C-type sesquioxide structure according 
to studies carried out by Zachariasen (1928) and by Pauling and Shappell 
(1930). 

The way in which Si** ions are accommodated in the Mn,O, structure 
in the mineral braunite has been the subject of a great deal of speculation and 
some investigations, In natural occurrences the composition of this mineral 
seems to be close to the formula 3Mn,0;-MnSiO;, corresponding to approxi- 
mately 10 Wt.% SiO., although slight variations in this ratio as well as pres- 
ence of other oxides in the structure are often encountered, Analytical data 
for this mineral are summarized in Dana’s System of Mineralogy (1944), The 
authors of this book are careful to point out that the exact mechanism of the 
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substitution is unknown. Goldschmidt (1954), however, states categorically 
that “braunite is not contaminated manganese sesquioxide but the compound 

Some recent magnetic and chemical studies have shed new light on this 
problem. Fyfe (1949) has developed a chemical test for the presence of 
trivalent manganese, based on the formation of a stable complex ion with 
acetylacetone. He found that a—Mn.O, gives the reaction of Mn?+ + Mn!‘ 
whereas y—Mn.Q,, studied for comparison, reacts as 2 Mn**. Krishnan and 
Banerjee (1939) have made magnetic measurements on natura! braunites to 
determine the oxidation state of manganese in this phase. Magnetic suscepti- 
bility measurements were non-conclusive, because the magnetic moment of 
Mn**, 4.90 Bohr magnetons, is almost midway between the moments of Mn** 
and Mn‘**, 5.92 and 3.87 Bohr magnetons, respectively. However, measure- 
ments of the magnetic anisotropy of the braunite revealed that the manganese 
was present as Mn + Mn‘** rather than as Mn**. Krishnan and Banerjee 
therefore suggest that the formula for braunite be written 3MnMnO,. MnSiO 
as against the usual formula 3Mn.0,-MnSiO3. 

The present investigation definitely establishes that a wide range of com- 
positions is possible in braunite, from pure Mn.O, as one end member to a 
silica content of at least 40 Wt.%o SiO. It appears likely that the mechanism 
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Fig. 2. Diagr howil ompositions of MneOs(ss) crystals in air at 1000°C 
(solid dots). The light straig lines is the join MnO-SiO:—Mn.O; (compare discussion 
in text) 
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for the solid solution formation is a substitution of Sit* for Mn** in the 
Mn.O, lattice. The present data give strong support for the assumption that 
Mn.O, consists of divalent and tetravalent manganese rather than irivalent 
manganese, for the continuous and extensive solid solution established for 
Mn.O,(ss) is most likely explained by assuming the formula Mn?+ Mn‘*+ Os. 

A further support of the assumption that the Mn.O;(ss) crystals can be 
considered a solid solution Mn?+(Mn‘*+Si*+)Os; is offered by the results of 
chemical analysis shown in figure 2. Compositions of MnO, (ss) crystals (solid 
dots) prepared by heating manganese oxide—SiO, mixtures in air for 200 
hours at 1000°C fall close to the join MnO-SiO.—Mn.0, as predicted by the 
mechanism of substitution suggested above. 

The crystalline phases tephroite (2MnO-SiO.), tridymite (SiO,), cris- 
tobalite (SiO,.) and Mn,O, appear to be pure phases corresponding to the 
chemical formulas written above, as closely as could be determined with the 
experimental techniques used (x-ray and microscopic examination as well as 
phase evidence). The phase designated MnO-SiO.(ss), after quenching to 
room temperature, is closely similar to the naturally occurring mineral rhodo- 
nite. A small excess of manganese over the stoichiometric 1:1 ratio is permitted 
in this structure, giving rise to a solid solution one-phase area (labeled 
MnO-SiO.(ss)) in the upper diagram in figure 1. No high temperature x-ray 
examination of this phase was carried out in the present investigation. In view 
of the complexity of the structure of metasilicates (compare for instance ex- 
tensive studies of the magnesium analogue, MgQ-SiO. (Atlas, 1952)), an 
x-ray study of the MnO-SiO.(ss) phase as a function of O. partial pressure 
and thermal history is a research problem in itself. In the present paper the 
phase is referred to as MnO-SiO.(ss), leaving open the question as to the 
exact crystal structure of this phase. 

Another notable feature of the phase relations in this system is the stable 
existence in air of tephroite as well as MnO-SiO.(ss) at liquidus temperatures, 
hoth phases containing essentially all manganese as Mn**. This situation is 
completely different from the one prevailing in the system iron oxide—SiQz, 
where magnetite (Fe,Q,) and silica (cristobalite or tridymite) are the only 
crystalline phases present in equilibrium at liquidus temperatures in air, In 
other words, Mn**+ is much more stable relative to the other valence states of 
manganese than Fe** is relative to Fe**+. The extreme stability of Mn** in 
oxide phases is of course inherent in the numerical values listed in the free 
energy tables compiled by Coughlin (1954). 

As temperature is lowered, a tendency for oxidation of Mn** to higher 


valence states prevails, because this is a strongly exothermic reaction, Tephro- 


ite oxidizes at temperatures only slightly below the lowest liquidus temperature 
in the system, first to Mn,Q, plus MnO-SiO.(ss), and subsequently to 
MnO, (ss). The MnO-SiO.(ss) phase persists in air down to a temperature of 
1048°C before oxidation to Mn.Qs(ss) and tridymite takes place. 

The boundary curve between the fields of tridymite and cristobalite has 
been dashed at 1470°C in accordance with the “classical” picture of stability 
relations among the various modifications (Fenner, 1913), and assuming no 
solid solution in tridymite or cristobalite. No consistent data regarding the 
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relative stabilities of these phases in the present system were obtained in runs 
of 2 days duration. Tridymite always formed at temperatures below 1480°C 
and cristobalite above this temperature in a mixture containing 40 Wt.% 
manganese oxide, while cristobalite was the phase always forming in the tem- 
perature interval 1430-1475°C in a mixture containing 8 Wt.% manganese 
oxide. 

The incongruent melting of tephroite in air as distinguished from the 
congruent melting found by Glasser (1958) for this phase under strongly re- 
ducing conditions is a noteworthy feature. It is probable that the discrepancies 
in previous literature regarding the melting relations of this phase are due to 
unsatisfactory methods for controlling and maintaining the desired level of O. 
pressures in the gas phase These relations are likely to he better understood 
upon completion of studies of the system MnO—Mn,.0,—SiO, at various O, 
partial pressures now in progress in our laboratories. 

The extremely strong stabilizing effect of SiO. on Mn*+ relative to the 
higher valence states of manganese in melts becomes apparent by an inspection 
of the composition curve for the liquid in the lower half of figure 1. This curve 
is for liquids at liquidus temperatures in air, hence temperature is a variable 
in addition to composition (silica to manganese oxide ratio). The effect of 
temperature is simple: Because the oxidation of Mn** is exothermic, a de- 
crease in temperature causes a decrease in the Mn?+/Mn‘** ratio, Referring 
to figure 1 it will be noticed that a pure manganese oxide liquid in air at tem- 
perature slightly above the liquidus consists of approximately 40 Wt.4¢ MnO 
and 60 Wt.“¢ Mn.O,. As SiO, is added, liquidus temperatures decrease (see 
upper half of fig. 1). If temperature were the only factor determining the 
Mn Mn equilibrium, this ratio would be expected to decrease as silica is 
added, corresponding to a curve originating at the 40 MnO 60 Wt.% Mn.O 
composition point and located on the Mn.O, side of the straight line connecting 
this point with the SiO, corner of the triangle in figure 1, The actual curve, 
however. bends sharply to the opposite side of this straight line. This situation 
is a result of the strong effect of acidity’ on the Mn*+/Mn** equilibrium. As 
the acidity of the system increases with increasing SiO. content. the Mn?+ is 
stabilized relative to Mn°~*, and so strong is this stabilizing effect of SiO. that 
it completely overshadows the opposite effect on the Mn**+ /Mn°* equilibrium 
of decreasing temperatures. Hence the 0.21 atm. O, isobar along the liquidus 
surface has the course represented by the heavy curve in the lower half of 
figure 1. The slope of this curve varies continuously, except where a boundary 
curve ts crossed 


The composition of the silica-rich liquid in the region of equilibrium 


coexistence of two immiscible liquids was only estimated roughly, The refrac- 
tive index of the glass was not significantly larger than that of pure SiO, glass, 
and the manganese oxide content is judged not to exceed 1 Wt.%. 

The relation between a diagram drawn with the appearance of a binary 
system, such as the join Mn,O0,—SiO, represented in the upper half of figure 
1. and the ternary system MnO—Mn.0,—SiO, represented in the lower, tri- 


For a discussion of acid-base relations in silicate liquids, see for instance a recent paper 


by Weyl (1956). 
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angular diagram in figure 1, has been explained in detail in previous papers 


by this author dealing with phase equilibria in iron oxide containing silicate 
systems (Muan, 1955, 1957, 1958). The same methods described in those 
papers have been used for projecting composition points from the lower ternary 
system up to the “binary” join in figure 1. Straight lines are drawn through 
composition points with a direction toward the O corner of the triangle repre- 
senting the system Mn—Si—O. These lines are “reaction lines” describing 
changes in total composition of condensed phases as these react with O, of the 
gas phase in response to temperature changes. The intersections between these 
reaction lines and the Mn,0,—SiO, join in the lower triangular diagram are 
taken to represent compositions of condensed phases in the upper simplified 
“binary” diagram, These points are projected to the upper, binary join by 
first swinging the Mn,O0,—SiO, join in the triangle around the SiO, corner 
upward until it falls on the Mn.O,—SiO, join, and then displacing the points 
along vertical lines to the upper half of the figure. 

When using and interpreting the “binary” phase diagram in the upper 
part of figure 1 it is important to keep in mind its relation to the ternary 
system. Although it is true that compositions of condensed phases cannot be 
expressed in terms of two chosen components, such as Mn,O, and SiQ., never- 
theless this diagram can be treated like a binary system as far as path of 
crystallization is concerned, The oxygen of the atmosphere (air) reacts with 
the condensed phases to keep the liquid compositions along the 0.21 atm. O, 
isobaric curve regardless of the nature of the crystalline phases separating out. 
Hence an apparently invariant situation will exist where the 0.21 atm, O, iso- 
bar crosses a boundary curve. Only at one temperature, at the given O, partial 
pressure, can two crystalline phases and a liquid coexist in equilibrium, be- 
cause one degree of freedom has been “used up” in choosing air as the at- 
mosphere. Consequently, intersections between the 0.21 atm, O, isobar and 
univariant lines in the system Mn—Si—O appear as invariant situations in the 
upper “binary” system. 

An extensive discussion of the application of phase equilibria in the 
system manganese oxide-silica to geochemical and technological problems 
must await results of studies presently in progress in our laboratories at 
various Q, partial pressures. However, a few general remarks may be made. 
The relatively constant ratio of silica to manganese oxide in naturally oc- 
curring braunites is an interesting observation which may perhaps be used to 
infer some of the conditions existing during the formation of this phase in 
nature. The composition limits of braunite depend primarily on two para- 
meters, temperature and Q, partial pressure. As O, partial pressure is reduced, 
the Mn.O, to Mn,O0, decomposition temperature decreases, and the two 
boundary curves in figure 1, representing compositions of braunite in equili- 
brium with Mn,O, and MnO-SiO.(ss), respectively, move downward, Hence 
the “peak” at which these lines meet also moves downward with decreasing O, 
pressure, and the composition of the braunite phase in terms of silica to 
manganese oxide ratio at this peak may also vary. These relations are illus- 
trated in the sketch in figure 3. If a naturally occurring braunite has been 
formed at a temperature close to the “peak” under the prevailing conditions of 
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Fig Sketch presented in order to indicate probable changes in composition range 
of MnsOs(ss) as Oz parti iressure is decreased below that of air (compare discussion 


In text) 


Q, partial pressure, its composition range should be very narrow. Moreover, if 
the shift in the “peak” composition is known as a function of temperature and 
Q. partial pressure from laboratory examinations, a close estimate of the values 
of these parameters at the time of formation of the mineral assemblage may 
be possible. Although slow reaction rates probably will prevent an experimental 
study of these reactions down to temperatures much below 1000°C, an attempt 
to map out the composition limits of the braunite phase volume seems like a 


worthwhile subject for future research. 


SUMMARY 
4 start has been made on the investigation of phase equilibria in the 
system Mn—Si—O by working along the 0.21 atm. O. isobaric section 
through this system at a total pressure of 1 atm. The phase relations are 
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represented in an apparently binary diagram (Mn,0,—SiO,) showing phases 
present as a function of temperature, and a supplementary diagram showing 
true compositions of condensed phases. Intersections between the 0.21 atm, O, 
isobar and boundary curves on the liquidus surface correspond to liquidus 
univariant situations in the ternary system characterized by the equilibrium 
coexistence of two crystalline, one liquid and a gas phase, or one crystalline, 
two liquid and a gas phase. However, because the O, partial pressure of the 
latter has been chosen, the existing degree of freedom has been “used up”, and 
the presence of these phases in equilibrium appears as invariant situations in 
the system drawn as the “binary” system Mn,0,—SiO, for sake of simplifica- 
tion. Similarly, subsolidus ternary univariant situations in the system Mn— 
Si—QO., with 3 crystalline phases and gas existing together in equilibrium, ap- 
pear as invariant situations in the simplified “binary” diagram. 

The “binary” system and the “invariant” situations are characterized as 
follows: manganese oxide melts at 1567°C in air, cubic Mn,0, being the pri- 
mary crystalline phase. Liquidus temperatures decrease sharply as SiO, is 
added A “peritectic” situation exists at 1230°C with Mn,0,, tephroite and a 
liquid of composition 54MnO, 13Mn.0,, 33 Wt.% SiO. present together in 
equilibrium with the gas. At 1206°C tephroite, MnO-SiO.(ss) (rhodonite), 
a liquid of composition 54MnO, 12Mn.0,, 34 Wt.% SiO, and gas coexist in 
equilibrium in what appears to be an eutectic situation. With SiO, contents in- 
creasing above 34 Wt.%, liquidus temperatures start rising, MnO-SiO,.(ss) 
being the primary crystalline phase in the composition range 34-42 Wt.% 
SiO.. At 1272°C another apparently peritectic sitaution prevails, with MnO- 
SiO. (ss), tridymite and liquid of composition 53MnO, 5Mn.O,, 42 Wt.% 
SiO. existing together in equilibrium in air. As SiO, content increases still 
farther, liquidus temperatures rise very steeply toward a plateau at 1700°C 
caused by the presence of two immiscible liquids over the composition range 
55-99 Wt.% SiO. at liquidus temperatures. The melting point of pure SiO, 
(cristobalite) is 1723°C. Phase transformation in Mn,O, and oxidation- 
reduction equilibria give rise to four apparently invariant situations in the 
subsolidus region of the system in the temperature range investigated. At 
1204°C Mn, 0,. tephroite and MnO-SiO.(ss) coexist in equilibrium in air, At 
1168°C the following crystalline phases exist together in equilibrium: cubic 
Mn.,O,, Mn.O,(ss) and MnO-SiO.(ss). At 1160°C another “invariant” situa- 
tion exists as the tetragonal and cubic forms of Mn,O, are present together 
with Mn.O,(ss) and gas. The phase MnO-SiO,.(ss) is stable in air down to 
1048°C. At this temperature it decomposes to Mn,O;(ss) and tridymite in an 
apparently eutectoid situation. 

The presence of tephroite and MnO-SiO.(ss) as the equilibrium phases 
at liquidus temperatures in air shows the extreme stability of the Mn*+ ion 
and accounts for its occurrence almost to the exclusion of manganese ions of 
higher valence states in igneous minerals and rocks. The Mn-O;(ss) crystals 
appearing as one of the equilibrium phases in the subsolidus temperature 
region seem to be the synthetic analogues of the naturally occurring mineral 
braunite, “3Mn.O,-MnSiO,”. A large range of manganese oxide to silica 
ratios in this phase has been found in the present investigation, strongly sug- 
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gesting a subsitution of Sit** for Mn** in the lattice of a—Mn.O,. In accord- 
ance with this observation and previous evidence from chemical as well as 
magnetic tests the e-form of Mn.Q, is probably best represented by the formu- 
la (Mn Mn**)O,, and the braunite formula should be written more ap- 


propriately Mn** {Mn**, Si‘* )O,. Compositions of natural braunites should 


depend mainly on temperature and QO, partial pressure existing during forma- 
tion of the mineral. Further laboratory work on composition limits for the 
Mn.O,(ss) phase as a function of temperature and QO, partial pressure, com- 
bined with chemical and petrographic analyses of braunite in manganese de- 
posits, may shed some new light on the conditions of formation of these 
deposits. 
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THE ANALCIME-JADEITE PHASE BOUNDARY: 
SOME INDIRECT DEDUCTIONS 
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ABSTRACT From a consideratior 


n of phase changes involving analcime some estimates 

have been made of the equilibrium pressures of the reaction: analcime jadeite + H.O. 
The results, although subject to large uncertainties, suggest that equilibrium pressures 

from existir erimental data are too high, particularly at the low temperature end of 

analcime confirm this conclusion and indicate 
at low temperatures. 


the systen 


Over the past lew years there has been a steady stream of papers con- 


cerned with the stability of mineral assemblages with total compositions cor- 
re sponding to and These compositions 
are involved in the reactions 
nepheline albite 2 jadeite 
albite jadeite + quartz. 

In a hydrous system the place of the zeolite analcime must be considered, 
and from phase relations involving this phase it is possible to make an estimate 
of the conditions of the reaction: 

analcime jadeite + water. 
At the present time calorimetric data exists for the phases albite, jadeite. 
nepheline and quartz but only entropy data is available for analcime. However. 


estimates of the necessar\ qu intities can be obtained from phase relations, 


The analcime nepheline albite water phase boundary and an estimate 
of the analcime jadeite equilibrium pressure at 535°C.—F rom studies in 
this laboratory (Valpy. 1958) it is considered that at 535°C (+ 10°) the 
equilibrium pressure of the reaction: 


analcime 1, nepheline + 1 albite + H.O (1) 
is LOOO atmospheres This value is in good agreement with that of Yoder 
(1954). At this point we then have the condition that the free energy change 
of reaction (1) is zero. From calorimetric and volume data the free energy of 
the reaction 

nepheling albite jadeite (2) 
be estimated at 535°C and 1000 atmospheres and hence under these con- 


ditions the free energy of reaction (2) is identical with the free energy of the 


may 


reaction 


inalcime jadeite + H.O (3) 
From the data summarized by Kelley. et al. (1955) the free energy of reaction 


(3) at 535°C and 1000 atmospheres is approximately 3,400 calories, ( All 


figures are rounded to the nearest 100 calories as errors involved in the figures 


themselves and approximations in the subsequent calculations do not justify 
greater accuracy.) As analcime is a phase of low density. at any moderate 
pressure the AV of reaction (3) is large and negative. In fact. if the density 


of water is unity. then the AV of reaction (3) is larger than that of reaction 


) 


16 
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To evaluate the equilibrium pressure of reaction (3) at 535°C we require 
data on the molar volume of water between 1000 atmospheres and this pres- 
sure. We have attempted to estimate this by extrapolation of Kennedy's data 
(1950) making use of the linear character of constant density, pressure- 
temperature plots. With this data the equilibrium pressure estimated is 10,700 
atmospheres. The average water density in the range was taken as 0.82 with an 
average AV of —14.43 em*, The density of analcime was taken as 2.267 from 
Hartwig (1930). Other densities were taken from Robertson, Birch and 
MacDonald (1957). Although the water density is difficult to estimate and may 
he considerably in error the AV due to this uncertainty should not be more 
than 15 percent and hence the estimated pressure should be within + 1500 
atmospheres. The AV of reaction (2) is —16.8 cm’, only slightly greater than 
AV of (3). and hence the transition pressure of (3) should not be much 
greater than that of (2). In these calculations thermal expansion and com- 
pressibility of all solid phases is ignored. 


The analcime + quartz albite + water phase boundary and an estimate of 
the analcime jadeite + water equilibrium pressure at 285°C.—The above 
reasoning may now be applied to the reaction: 
analcime + quartz albite + H.O (4) 
which we consider to be in equilibrium at 285°C (+ 15°) and 1000 atmos- 
pheres (Ellis, 1958; Valpy, 1958). From Kelley. et al. (1953) thermodynamic 
data is available for the reaction: 
albite jadeite + quartz (5) 
and as the free energy change of reaction (4) is zero at 285° and 1000 atmos- 
pheres the free energy change of reaction (5) will be identical under these 
conditions with that of the reaction: 
analeime + quartz jadeite + quartz + H.O (6) 
The free energy change of (6) is thus approximately 4.200 calories, This free 
energy is greater than that of reaction (3) at 535° and 1000 atmospheres but 
at 285°C the AV will be greater due to the smaller molar volume of water. In 
this region extrapolation of Kennedy's data was subject to considerable un- 
certainties and an average water density of unity was assumed, This figure is 
possibly a little high but should not introduce more than 10 percent error in 
AV. The average AV of reaction (6) used was —18.4 cm‘ and the equilibrium 


pressure of reaction (6) is thus 10,500 atmospheres, The pressure is very 


similar to that estimated at 535°C and indicates that reaction (3) (or (6), 
which is of course, identical) has a very small entropy change. 


The equilibrium pressure of the reaction: analcime jadeite + water at 
To obtain an estimate of the transition pressure of analcime to jadeite 

C it is necessary to have data on the entropy of analcime. This is avail- 

able from the measurements of King (1955) who gives a value of 56.0 cals 
mol./degree at 25°C. This is an abnormally high entropy. Following Latimer 
(1951) it might be expected that the entropy would be close to that given by 
1, San + % + 19.5 cals./mol./degree. In this calculation 
Si. is taken as 9.6 cals. mol.. a value obtained from cousideration of a large 
number of hydrates. The entropy is thus 6.5 cals./mol. too large, In general, 
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there is a good correlation between entropy and molar volume (see Fyfe, 
lurner and Verhoogen, 1958) and relative to nepheline, albite and ice, the 
molar volume of analcime is similar. This suggests that the large entropy may 
be associated with randomness of the cation arrangement in the analcime 
structure. The entropy and molar volume of jadeite are lower than that one 
would estimate by summing values for nepheline and albite. It may be noted 
further that there is a general correlation between the entropy contribution of 
water to a solid and the apparent volume of water in the solid, (this is taken 
as the difference in the volume of hydrate and anhydrous phase). In many 
cases the apparent molar volume of water is in the range 13-17 cm* with 


entropy contributions ranging from 8-11 cals./mol. With analcime the apparent 


molar volume of water is 20.3 cm‘ and a high entropy contribution would be 
expected. A further confirmation of the high entropy of analcime can be ob- 
tained from the temperature interval between the curves of reactions (4) and 

At 535°C and a thousand bars the free energy of reaction (4) is identical 
with that of the reaction: 


nepheline quartz 1, albite 
at 535°C. The free energy of (7) is available from Kelley, et al. (1953). If we 
assume that AC, for reaction (4) is small, then AT AS AG where AT is the 
difference in boundary temperatures of reactions (1) and (4), 250°C, and AS 
is the average entropy of reaction (4) between 285° and 535°C. From this 
argument AS is 11.6 cals./mol. This is a low value for a dehydration reaction 
in the region 285-535°C (see Fyfe, Turner and Verhoogen, 1958). From data 
in Kelley (1949) and from estimates of the entropy of water by Fyfe, Turner 
and Verhoogen, we may compare this value with those for the reactions brucite 
to periclase + vapor and boehmite to corundum + vapor. The average AS 
for these reactions is near 14 cals./mol. Using 56.0 cals./mol. for the entropy 
of analcime the entropy change of reaction (3) at 25°C is —7.35 cals./mol. 
degree. The volume change is certainly negative so that the low temperature 


200 300 400 500 600 


Fig. 1. Some estimated boundary relations concerned with jadeite, F-B and A-E are 
the calculated curves for reactions (2) and (5). A-B is the estimated curve for reaction 
(3). D-C is the boundary for (3) reported by Griggs and Kennedy. H and G represent 
the experimental boundaries for reactions (1) and (4) from Ellis, 1958, and Valpy, 1958. 
Dotted lines represent approximate corrections from data of Robertson, Birch and Mac- 
Donald (1957) 
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boundary should have positive slope. This conclusion is not changed by con- 
sideration of the compressibility of analcime (Yoder, 1955). 

We may now proceed to estimate the free energy change of reaction (4) 
at 25°C. A “normal” AS of hydration in the range 25-285°C would be around 
8 cals./mol. (Fyfe, Turner, Verhoogen, 1958) and for analcime it is more like- 
ly to be around 2-3 cals./mol. The entropy of reaction (4) at 25°C has the 
remarkably small value of + 1cal./mol. These values lead to a free energy 
change at 25°C and 1000 bars for reaction (4) of 780 cals./mol. and for re- 
action (3) a free energy change of 2,900 cals./mol, Assuming a water density 
of unity the equilibrium pressure would be 7,600 atmospheres and allowing for 
the compression of water (Dorsey, 1940), is more likely to be around 7000 
atmospheres. In figure 1, data for all calculated boundaries are summarized. 
From the AS and AV values the positive slope of the analcime-jadeite boundary 
should be around 15 atmospheres/°C at 25°C, 

It should perhaps be mentioned that a similar value at 25°C is obtained 
simply by taking an average entropy of reaction (1) from 535-25°C assuming 
that the entropy of analcime at 25°C is 56 cals./mol./degree. This then frees 
the calculations from any use of data from reaction (5). 

DISCUSSION 

The only experimental boundary between the phases analcime and jadeite 
is that in Griggs and Kennedy (1956). This boundary is shown in figure 1 line 
C-D. Obviously there is enormous departure from the present estimated curve 
A-B in the low temperature region. Although there are great uncertainties in 
the estimated curve it is difficult to account for a three-fold error. The experi- 
mental results of Robertson, Birch and MacDonald (1957) on jadeite stability 
suggest that the calculated boundaries A-E and B-F are about 2000 bars too 
low and if these are raised the analcime-jadeite boundary will be increased by 
a similar amount. This would lead to close agreement with the published 
boundary at 600°C, and the position of the invariant point (B’) indicated by 
Robertson, Birch and MacDonald. 

It will be noticed that the position of the second invariant point (A) sug- 
gests a large amount of back bending of the analcime-albite boundary. Such 
back bending must occur in this boundary because AV of the reaction will 
change sign at high pressures but it is difficult to estimate the extent. On the 
other hand, if line C-D was correct then the invariant point would be at C and 
there would be little scope for the required reversal of slope which once more 
suggests that this boundary is too high. It is also possible that the observed 
analcime-albite boundary is too high, it certainly represents an upper limit of 
analcime rather than a lower limit of albite as will be discussed in a later paper. 

Features of petrogenetic interest arise from the predicted flat analcime- 
jadeite boundary. There is some evidence that jadeite is formed in a moderate- 
ly low temperature environment. (e.g. Yoder and Chesterman, 1951; Coleman, 
1954), with hydrous minerals. The formation of low temperature jadeite pre- 
sents no problems if the total pressure in the environment is much greater than 
the partial pressure of water, but in the glaucophane schists of California, 
jadeite is associated with lawsonite, a mineral which, judging from present 
data of one of the authors (W.S.F.), requires a considerable partial pressure 
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of water to form. (A limited number of experiments indicate pressures greater 
than 5000 bars at 400°C.) Also lawsonite is unusual in that its molar volume 
is similar to that of anorthite so that differential pressure has little influence on 
the lawsonite-feldspar equilibrium. Thus high water pressures may exist in the 
region where jadeite forms in association with lawsonite. If the Griggs-Kennedy 
boundary was correct, the depth of burial required in a low temperature en- 
vironment, where water and load pressures are equal, would be so great as to 
suggest that jadeite could not form in this region. But the predicted curve 
would possibly be attained in regions of lowest thermal gradients (see Birch, 
1955). where at 30 kilometers the temperature might be less than 300°C. 

All these remarks require experimental confirmation, particularly in view 


of possible effects of solid solutions. In general these will lower pressures neces- 


sary for the formation of impure jadeite. 
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